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This  dissertation  mainly  focuses  on  analysis  and  modeling  of 
parasitic  effects  in  advanced  silicon-on-insulator  (SOI)  complementary 
metal-oxide-semiconductor  (CMOS)  technology.  As  the  channel  length  is 
reduced  to  the  nanometer  regime,  the  suppression  of  parasitic  effects 
becomes  an  increasingly  difficult  technological  challenge.  First,  the  direct 
gate  tunneling  current  is  a major  parasitic  leakage  component  because 
extremely  thin  gate  oxide  is  required  in  order  to  provide  sufficient  current 
drive  with  reduced  supply  voltage,  and  to  control  short-channel  effects 
(SCEs).  Gate-to-body  tunneling  current,  which  modulates  the  body  voltage 
in  scaled  partially  depleted  (PD)  SOI  MOS  field-effect  transistors  (FETs), 
is  physically  modeled  and  implemented  in  the  University  of  Florida  PD 
SOI/Bulk  (UFPDB)  MOSFET  model.  Several  scaling  effects  (i.e., 
quantization  effect  in  the  channel,  exchange  energy  of  inversion  electrons, 
and  Fermi-Dirac  statistics  in  the  high-doped  gate)  are  identified  and 


accounted  for.  Using  the  model  in  UFPDB,  it  is  shown  how  the  tunneling 
current  tends  to  suppress  dynamic  floating-body  effects  in  some  PD  SOI 
circuits,  providing  benefits  that  follow  from  the  complex  interactions 
among  the  bias  conditions,  circuit  topologies,  and  switching  patterns.  The 
parasitic  bipolar-junction-transistor  (BJT)  effects  in  the  floating-body  SOI 
MOSFET  are  modeled  because  they  can  induce  transient  BJT  current  and 
amplify  the  gate-induced-drain-leakage  (GIDL)  current  in  nonclassical 
devices.  The  previous  version  of  the  BJT  model  in  the  University  of  Florida 
double-gate  (UFDG)  MOSFET  model  is  modified  and  upgraded  to  describe 
correctly  the  majority-carrier  concentration,  which  governs  the  parasitic 
BJT  gain,  in  the  (undoped)  ultra-thin  body  (UTB)  of  nonclassical 
MOSFETs.  The  triple-gate  MOSFET  (TGFET),  proposed  as  a means  to 
alleviate  the  stringent  UTB  processing  requirements,  can  be  undermined 
by  the  top  corners  of  the  silicon-fin  body,  which  can  significantly  affect  the 
channel  current-voltage  characteristics  of  the  TGFET.  It  is  shown  that  the 
corner  effects  can  be  eliminated  by  leaving  the  body  undoped.  Then,  for 
viable  undoped  bodies,  the  UTB  dimensions  needed  to  control  SCEs  in 
nanoscale  TGFETs  are  examined  by  three-dimensional  numerical  device 
simulations.  The  results  show  that  much  more  stringent  body  scaling  is 
needed  for  undoped  TGFETs  relative  to  the  doped  ones,  which  are 
technologically  and  electrically  infeasible.  When  the  undoped-body 
dimensions  are  scaled  for  adequate  SCE  control,  it  is  found  that  the 
TGFET  suffers  from  a significant  layout-area  disadvantage  relative  to  the 
double-gate  (DG)  FinFET  and  the  planar  fully-depleted  MOSFET 
(FDFET),  and  thus  it  is  not  feasible.  It  is  concluded  that  the  DG  FinFET, 
which  is  more  scalable  than  the  FDFET,  is  the  nonclassical  device  with  the 
most  potential  for  future  nanoscale  CMOS  applications. 


IX 


CHAPTER  1 
INTRODUCTION 


Scale-down  of  complementary  metal-oxide-semiconductor  (CMOS) 
device  dimensions  in  conventional  bulk-silicon  CMOS  technology  has 
been  a primary  driver  over  the  past  three  decades  in  order  to  attain 
continued  improvement  in  integrated-circuit  (IC)  performance  and 
reduction  in  size  (for  lower  cost).  Silicon-on-insulator  (SOI)  technology 
has  emerged  as  a competitor  with  bulk-silicon  technology  for  high 
performance  and  as  an  effective  means  of  extending  CMOS  scaling 
beyond  bulk-silicon  technology  limits.  Now,  the  partially  depleted  (PD) 
SOI  MOS-field-effect-transistor  (MOSFET),  which  is  already  in  product 
stage,  is  rapidly  replacing  the  bulk-silicon  MOSFET  for  high- 
performance  applications;  and,  with  ultra-thin-bodies  (UTB),  the  fully 
depleted  SOI  MOSFET  and  multi-gate  MOSFETs  are  most  promising 
candidates  for  scaling  beyond  the  limit  of  conventional  bulk-silicon 
CMOS  technology. 

As  the  channel  length  is  reduced  to  the  nanometer  regime, 
however,  the  suppression  of  parasitic  effects  becomes  an  increasingly 
difficult  technological  challenge.  Figure  1.1  shows  the  several  parasitic 
effects  of  classical  and  nonclassical  CMOS  device  structures.  First,  the 
direct  gate-tunneling  current  can  be  a major  parasitic  leakage  component 
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Figure  1.1  The  several  parasitic  effects  of  classical  and  nonclassical 
CMOS  device  structures,  with  the  implied  works  done  in  this 
dissertation  indicated.  Also,  the  structures  of  the  primary 
nonclassical  CMOS  devices  are  illustrated. 
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because  extremely  thin  gate  oxide  is  required  in  order  to  provide 
sufficient  current  drive  with  reduced  supply  voltage,  and  to  control  short- 
channel  effects  (SCEs).  Of  the  several  components  of  the  direct  gate- 
tunneling current,  the  gate-to-body  tunneling  current  modulates  the  body 
voltage  and  can  influence  the  floating-body  effects  in  scaled  PD  SOI 
MOSFETs.  Therefore,  the  modeling  of  the  gate-to-body  direct  tunneling 
current  and  assessment  of  its  effects  are  crucial  for  the  design  of  floating- 
body  (FB)  PD  SOI  circuits.  The  parasitic  bipolar-junction-transistor 
(BJT)  effects  of  the  FB  SOI  MOSFET  are  also  important  because  they 
induce  transient  BJT  current  and  amplify  the  gate-induced-drain- 
leakage  (GIDL)  current.  Thus,  modeling  of  the  parasitic  BJT  in  undoped 
UTB  devices  illustrated  in  Figure  1.1,  the  most  viable  ones  because  of 
technological  doping  issues,  is  required  to  simulate  precisely  the  dynamic 
and  static  off-state  leakage  currents  in  view  of  transient  circuit  operation 
and  power  consumption.  The  triple-gate  MOSFET  (TGFET)  illustrated  in 
Figure  1.1,  which  is  proposed  as  a means  to  alleviate  the  stringent  thin 
silicon-fin  requirements  of  UTB  devices,  is  subject  to  excessive 
conduction  in  the  top  corners  of  the  silicon-fin,  which  can  significantly 
affect  the  channel  current-voltage  characteristics  of  the  TGFET.  Analysis 
of  the  parasitic  corner  effects  is  needed  to  examine  the  feasibility  of  the 
TGFET. 

As  CMOS  technology  is  scaled  to  the  nanometer  gate-length  regime 
for  higher  density  and  speed,  thinner  gate  oxide  is  required  in  order  to 
provide  sufficient  current  drive  with  reduced  supply  voltage,  and  to 
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control  short-channel  effects.  Sub-lOOnm  CMOS  technologies  will  require 
the  (equivalent)  gate  oxide  thickness  to  approach  ~lnm  [SemOl].  As  the 
actual  thickness  (tox)  is  decreased,  the  direct  gate-tunneling  current 
increases  with  nearly  exponential  dependence  on  it.  The  higher  gate 
current  not  only  increases  the  standby  power  consumption  of  a highly 
integrated  CMOS  circuit,  but  also  can  adversely  affect  MOSFET 
performance.  Therefore,  previous  studies  of  direct  tunneling  currents  in 
MOS  devices  have  focused  primarily  on  the  predominant  component 
which  flows  between  the  gate  and  the  channel  [Reg99,  Yan99,  KhaOO, 
Lo97,  YanOO];  the  much  smaller  tunneling  current  that  flows  between  the 
gate  and  the  body  has  not  been  investigated  extensively  [MajOO,  Sha99, 
LeeOO],  nor  modeled  properly.  However,  in  PD  SOI  MOSFETs,  the  latter 
current  modulates  the  body  voltage  and  can  influence  the  floating-body 
effects  [FunOO,  JosOl,  Chu02,  Poi02],  Such  influence,  which  can  be 
beneficial  or  detrimental,  can  prevail  even  for  relatively  thick  values  of 
tox  for  which  the  larger  gate-channel  current  is  negligible.  In  chapter  2, 
we  analyze  the  gate-to-body  direct  tunneling  current  based  initially  on 
Harrison’s  independent-electron  tunneling  model  [Har61].  However,  for 
the  inversion  case,  the  proper  characterization  of  the  gate-to-body 
tunneling  current  is  somewhat  different  from  Harrison’s  model  since  the 
electron  density-of-states  effective  mass  is  different  in  the  valence  and 
conduction  bands,  which  contain  the  initial  and  final  states  of  the 
tunneling  electrons  that  constitute  the  main  component  of  the  current. 
Moreover,  several  scaling  effects  (i.e. , quantization  in  the  channel, 
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exchange  energy  of  inversion  electrons,  and  Fermi-Dirac  carrier 
distribution  in  the  poly-Si  gate)  and  impact  ionization,  which  have  been 
overlooked  in  contemporary  modeling  of  gate  tunneling  current,  are 
identified  and  accounted  for.  We  generalize  the  presentation  by 
discussing  the  several  gate-to-body  tunneling  current  components  in  both 
the  nMOSFET  and  the  pMOSFET,  for  both  inversion  and  accumulation 
cases,  after  which  we  describe  our  physics-based  analytic  modeling  of  the 
predominant  components.  The  model  is  implemented  in  our  process/ 
physics-based  compact  model  (UFPDB  [Fos02a]),  which  is  unified  for  PD 
SOI  and  bulk-silicon  MOSFETs. 

The  FB  effects  in  PD  SOI  MOSFETs,  and  the  resulting  hysteresis 
and  instability  in  dynamic  operations,  pose  major  challenges  for  very- 
large-scale-integration  (VLSI)  circuit  designs.  Further,  the  gate-to-body 
tunneling  current  modulates  the  body  voltage,  which  is  varied  by 
capacitive  coupling  in  transient  operation,  and  thus  can  influence  the  FB 
effects  [FunOO,  JosOl,  Chu02,  Poi02].  Such  influence  can  be  beneficial  or 
detrimental  for  the  dynamic  operations  of  FB  PD  SOI  circuits  and  must 
be  rigorously  examined  by  a reliable  physical  model.  In  chapter  3,  using 
the  gate-to-body  current  model  in  UFPDB/Spice3,  we  assess  the 
significance  of  the  gate-to-body  current  for  circuit  topologies  and 
switching  patterns  susceptible  to  FB  effects. 

The  BJT  effects  of  the  FB  SOI  MOSFET  have  constituted  a critical 
issue  for  a long  time  because  they  induce  a premature  breakdown  [You88, 
Cho91],  latch  [Che88,  Cho91],  transient  BJT  current  [Pel95],  and  also 
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amplify  the  GIDL  current  [Che92].  While  impact  ionization-induced 
phenomena  (premature  breakdown  and  latch)  have  become  mitigated  due 
to  scaled  operation  voltage,  the  transient  BJT  current  and  the  amplified 
GIDL  current  are  still  major  concerns  for  nonclassical  UTB  devices 
regarding  control  of  off-state  leakage  current  in  view  of  dynamic  and 
static  power  consumption  [Sak03,  Fos99,  Sch04].  However,  the  parasitic 
BJT  model  [Kri96]  of  UFDG  [Fos03a],  which  is  our  process/physics-based 
compact  model  for  UTB  devices,  is  inadequate  to  describe  the  transient 
BJT  action  and  the  amplification  of  GIDL  current  in  nonclassical  UTB 
devices  because  it  does  not  correctly  account  for  accumulation  in  the 
undoped  UTBs.  In  chapter  4,  the  previous  BJT  model  is  modified  and 
upgraded  to  correctly  describe  the  majority-carrier  concentration  in  the 
UTB,  which  governs  the  parasitic  BJT  gain.  Further,  circuit  application 
with  the  upgraded  UFDG  in  Spice3  [Fos03a]  is  exemplified  by 
simulations  of  a nonclassical  dynamic  logic  circuit,  which  show  the 
parasitic  BJT  effect  on  circuit  operation  and  stress  the  need  to  account 
for  it  in  IC  design  at  future  CMOS  technology  nodes. 

Conventional  scaling  of  classical  bulk-silicon  and  PD  SOI  CMOS 
has  become  very  challenging  because  control  of  SCEs  for  acceptable  I0i/Ioff 
ratios  requires  precise  channel  doping  levels  and  gradients  [FraOl]  that 
are,  ultimately,  impossible  to  achieve.  As  noted,  promising  alternatives 
for  continued  CMOS  scaling  are  the  nonclassical  devices  having  UTBs, 
i.e.,  the  planar  fully  depleted  single-gate  SOI  MOSFET  (FDFET)  [Dor02] 
and  the  quasi-planar  double-gate  (DG)  FinFET  [Ked02],  which  give  good 
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SCE  control  via  the  UTB.  However,  the  good  SCE  control  in  the 
nonclassical  nanoscale  devices  requires  that  the  UTB  thickness  (tgj)  be 
considerably  thinner  than  the  gate  length  (Lg),  as  well  as  being  uniform 
across  the  device  and  circuit.  The  TGFET  was  recently  proposed  [Doy03a, 
Doy03b]  as  a means  to  alleviate  the  stringent  thin-tSi  requirements  of 
UTB  FDFETs  and  DG  FinFETs.  However,  design  studies  of  TGFETs 
[Doy03b]  have  revealed  that  the  top  corners  of  the  silicon  body  can,  with 
strong  dependence  on  their  shape,  or  radius  of  curvature,  significantly 
affect  the  channel  current-voltage  characteristics  of  the  TGFETs. 
Typically,  the  corner  components  of  current  reflect  a lower  threshold 
voltage  (Vt)  and  higher  off-state  current  (I0ff)  relative  to  the  respective 
properties  of  the  bulk  TG  device,  and  for  some  body  structures  they  can 
cause  a “hump”  in  the  subthreshold  log(IDs)-VGS  characteristic  [Doy03b]. 
Hence,  a serious  technological  issue  is  implied.  And,  the  layout  area  of  the 
TGFET  is  an  issue  for  nanoscale  CMOS  [Doy03b].  The  feasibility  of  the 
TGFET  is  thus  not  yet  clear.  In  chapter  5,  we  examine  the  abnormal 
corner  effects  in  nanoscale  TGFETs,  using  two-dimensional  (2-D) 
numerical  device  simulations  and  quasi-2-D  analysis  to  gain  physical 
insights  on  how  the  effects  could  possibly  be  suppressed,  irrespective  of 
the  radius  of  curvature,  which  is  difficult  to  control  technologically.  Then, 
we  examine  the  feasibility  of  the  TGFET  for  nanoscale  CMOS  applications 
with  regard  to  SCEs  and  layout  area.  We  check  the  SCEs  for  varying 
dimensions  of  the  body,  with  and  without  doping,  using  a three- 
dimensional  (3-D)  numerical  device  simulator.  Then,  with  the  SCEs 
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adequately  controlled  via  body  scaling  in  accord  with  the  3-D  simulation 
results,  we  analyze  the  gate  layout  area  of  the  integrated  TGFET  needed 
for  current  drive,  and  compare  it  with  those  of  the  FDFET  and  DG 
FinFET. 


CHAPTER  2 

GATE-TO-BODY  DIRECT  TUNNELING  CURRENT:  MODELING  AND 

IMPLEMENTATION 

2.1  Introduction 

As  CMOS  technology  is  scaled  to  the  nanometer  gate-length 
regime  for  higher  density  and  speed,  thinner  gate  oxide  is  required  in 
order  to  provide  sufficient  current  drive  with  reduced  supply  voltage,  and 
to  control  SCEs.  Sub-lOOnm  CMOS  technologies  will  require  the 
(equivalent)  gate  oxide  thickness  to  approach  -lnm  [SemOl].  As  the  actual 
thickness  (tox)  is  decreased,  the  direct  gate-tunneling  current  increases 
with  nearly  exponential  dependence  on  it.  The  higher  gate  current  not 
only  increases  the  standby  power  consumption  of  a highly  integrated 
CMOS  circuit,  but  also  can  adversely  affect  MOSFET  performance. 
Therefore,  previous  studies  of  direct  tunneling  currents  in  MOS  devices 
have  focused  primarily  on  the  predominant  component  which  flows 
between  the  gate  and  the  channel  [Reg99,  Yan99,  KhaOO,  Lo97,  YanOO]; 
the  much  smaller  tunneling  current  that  flows  between  the  gate  and  the 
body  has  not  been  investigated  extensively  [MajOO,  Sha99,  LeeOO],  nor 
modeled  properly.  However,  in  PD  SOI  MOSFETs,  the  latter  current 
modulates  the  body  voltage  and  can  influence  the  floating-body  effects 
[FunOO,  JosOl,  Chu02,  Poi02].  Such  influence,  which  can  be  beneficial  or 
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detrimental,  can  prevail  even  for  relatively  thick  values  of  tox  for  which 
the  larger  gate-channel  current  is  negligible. 

In  this  chapter,  we  analyze  the  (poly-Si)  gate-to-body  direct 
tunneling  current  based  initially  on  Harrison’s  independent-electron 
tunneling  model  [Har61].  However,  for  the  inversion  case,  the  proper 
characterization  of  the  gate-to-body  tunneling  current  is  somewhat 
different  from  Harrison’s  model  since  the  electron  density-of-states 
effective  mass  is  different  in  the  valence  and  conduction  bands,  which 
contain  the  initial  and  final  states  of  the  tunneling  electrons  that 
constitute  the  main  component  of  the  current.  Moreover,  several  scaling 
effects  (i.e.,  quantization  in  the  channel,  exchange  energy  of  inversion 
electrons,  and  Fermi-Dirac  carrier  distribution  in  the  poly-Si  gate)  which 
have  been  overlooked  in  contemporary  modeling  of  gate  tunneling  current 
are  identified  and  accounted  for.  We  generalize  the  presentation  by 
discussing  the  several  gate-to-body  tunneling  current  components  in  both 
the  nMOSFET  and  the  pMOSFET,  for  both  inversion  and  accumulation 
cases,  after  which  we  describe  our  physics-based  analytic  modeling  of  the 
predominant  components.  The  model  is  implemented  in  our  process-based 
compact  model  (UFPDB  [Fos02a]),  which  is  unified  for  PD  SOI  and  bulk- 
silicon  MOSFETs. 

2.2  Direct  Gate  Tunneling  Current  Components 
The  components  of  gate  tunneling  current  are  identified  using 
the  simple  energy-band  diagrams  of  the  n+poly-Si/Si02/p-Si  nMOSFET 
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and  p+poly-Si/Si02/n-Si  pMOSFET  structures  in  Figure  2.1  and  Figure 
2.2,  respectively.  For  the  nMOSFET  biased  to  inversion  (VGB  > Vt)  in 
Figure  2.1(a),  the  conduction-band  electron  tunneling  ((D)  is  the  largest 
component  because  the  associated  Si-Si02  barrier  height  (3.1eV)  is  small 
compared  to  those  (4.2eV  and  4.5eV,  respectively)  associated  with  the 
valence-band  electron  tunneling  ((D)  and  the  valence-band  hole  tunneling 
((D).  However,  even  though  the  valence-band  tunneling  components  are 
relatively  small,  they  can  significantly  influence  FB  effects  in  PD  SOI 
nMOSFETs  because  they  result  in  hole  injection  in  the  body.  The  latter 
component  (D  is  seemingly  very  small  because  of  the  negligible 
concentration  of  holes  in  the  valence  band  of  the  gate.  However,  as  we 
suggest  herein,  it  becomes  predominant  at  high  VGB  because  of  the  impact 
ionization  (illustrated  in  the  figure)  driven  by  hot  electrons  in  the  gate 
from  component  (D,  which  enhances  the  hole  density.  Therefore,  all  three 
tunneling  components  in  Figure  2.1(a)  can  be  important  in  this  case. 

In  transient  operation  of  the  SOI  nMOSFET,  the  dynamic  VGB 
can  be  negative  due  to  the  capacitive  coupling  between  the  body  and  the 
source/drain.  Therefore,  the  accumulation  case  in  Figure  2.1(b)  must  be 
considered.  In  this  case,  all  three  gate  tunneling  current  components 
contribute  to  floating-body  effects  since  the  source  and  drain  are  isolated. 
The  conduction-band  electron  tunneling  (®)  is  the  dominant  component 
since  the  associated  Si-Si02  barrier  height  is  relatively  small.  The  impact 
ionization  driven  by  hot  electrons  in  the  substrate  and  hot  holes  in  the 
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Figure  2.1  Energy  bands  and  components  of  direct  gate  tunneling 

current  in  the  nMOSFET  with  n+  poly-Si  gate  for  the 
inversion  and  accumulation  cases.  The  impact-ionization 
effect  is  indicated  for  the  inversion  case.  Epg  and  Ep(j 
represent  the  Fermi-level  energy  (Ep)  of  the  silicon  and  the 
gate,  respectively. 
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Figure  2.2  Energy  bands  and  components  of  direct  gate  tunneling 
current  in  the  pMOSFET  with  p+  poly-Si  gate  for  the 
inversion  and  accumulation  cases.  The  impact-ionization 
effect  is  indicated  for  the  inversion  case. 
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gate,  from  components  ® and  ©,  respectively,  is  not  significant  because 
the  hole  and  electron  densities,  respectively,  in  the  substrate  valence  band 
and  the  gate  conduction  band  are  high  for  this  bias  condition.  Therefore, 
the  tunneling  component  ® is  the  only  important  one  in  this  case. 

For  the  pMOSFET  biased  to  inversion,  all  tunneling  components 
in  Figure  2.2(a)  are  analogous  to  those  of  the  nMOSFET  in  Figure  2.1(a), 
and  can  be  important.  However,  there  is  a difference  in  the  accumulation 
case  in  Figure  2.2(b).  Because  of  the  small  Si-Si02  barrier  height  of  the 
conduction-band  electron  tunneling  (©)  relative  to  the  valence-band  hole 
tunneling  (®),  the  component  © becomes  important  for  VGB  > VFB  when 
the  electron  density  in  the  substrate  is  substantial.  Therefore,  the  gate-to- 
body  direct  tunneling  current  is  mainly  the  component  ©,  whereas  it  is  the 
component  ® for  the  nMOSFET  in  Figure  2.1(b). 

In  summary,  the  gate-to-body  direct  tunneling  currents  are 
mainly  due  to  IVBE  (valence-band  electron  tunneling  © ) and  IyBH  (valence- 
band  hole  tunneling  ®)  at  inversion,  and  ICBE  (conduction-band  electron 
tunneling  ®)  at  accumulation  in  the  nMOSFET,  and  due  to  IyBE  (valence- 
band  electron  tunneling  ®)  and  ICBE  (conduction-band  electron  tunneling 
®)  at  inversion,  and  ICBE  (conduction-band  electron  tunneling  ©)  at 
accumulation  in  the  pMOSFET.  These  tunneling  currents  modulate  the 
threshold  voltage  by  changing  the  body  voltage  in  PD  SOI  MOSFETs,  and 
therefore  influence  the  FB  effects  [FunOO].  Herein  we  describe  physics- 
based  analytic  modeling  of  the  gate-to-body  direct  tunneling  currents,  and 
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its  implementation  in  our  compact  process-based  model  (UFPDB  [Pel02]) 
for  PD  SOI  and  bulk-silicon  MOSFETs. 

2.3  Inversion-Region  Model 
2.3.1  Valence-Band  Electron  Tunneling 

We  focus  first  on  the  direct  valence-band  electron  tunneling, 
which  is  a major  component  of  the  gate-to-body  current  at  the  normal 
inversion  mode  of  both  the  nMOSFET  and  the  pMOSFET  (before  the  onset 
of  impact  ionization  due  to  the  tunneling  from  the  channel).  As  seen  in 
Figure  2.1(a)  and  Figure  2.2(a),  the  mechanisms  are  exactly  the  same  in 
both  devices.  Therefore,  the  model  we  develop  here  for  the  nMOSFET  is 
directly  applicable  to  the  pMOSFET  as  well.  The  valence-band  electron 
tunneling  in  the  nMOSFET  for  VGB  > 0,  i.e.,  component  © in  Figure  2.1, 
as  governed  by  the  electric  potential  band  bending  \j/sb  in  the  silicon 
substrate  and  the  potential  drop  V|/p  in  the  gate  due  to  poly-Si  depletion 
[Tau98].  Valence-band  electrons  in  the  substrate  tunnel  into  the 
conduction  band  of  the  gate  when  EGg  (Eq  at  the  poly-Si  surface)  in  the 
gate  drops  below  Evs  (Ey  at  the  silicon  surface)  in  the  substrate,  as 
depicted  in  Figure  2.3.  Note  that  this  bias  condition  implies  strong 
inversion.  The  tunneling  current  density  J (in  x)  is  characterized  by 
modifying  the  independent-electron  tunneling  model  of  Harrison  [Har61] 
to  account  for  the  different  transverse  density-of-states  (DOS)  effective 
masses  in  the  silicon  valence  band  (mta*)  and  in  the  poly-Si  conduction 
band  (mtb*).  The  more  general  expression  for  the  net  (silicon  valence  band 
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(a)  Vgb  = 0 


(b)  VGB  - V0nset  <EVS  ~ ECg)  (c)  VGB  > VonSet  (EVS  " ECG  > 


Figure  2.3  Energy  bands  showing  the  onset  of  the  valence-band  electron 
tunneling  in  the  nMOSFET  (with  VBg  = 0). 
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(a)  to  poly-Si  conduction  band  (b))  electron  tunneling  current  density 
(magnitude)  is  given  by  (see  Appendix) 

Ja^b  = ^[mt<sfo"---'foPfaa-fb)dExdEt 

Pfbil  -f°)dE‘dE‘] ;(21) 

mta*,  mtb*:  transverse  DOS  effective  masses  in  the  states  a and  b, 
respectively; 

Ex:  energy  (kinetic  plus  potential)  in  the  x-direction  (normal  to  the  Si-Si02 
interface); 

Et:  transverse  (parallel  to  the  Si-Si02  interface)  energy  (kinetic) 
component; 

fa,  fb:  probabilities  (Fermi-Dirac)  of  occupation  of  states  a and  b, 
respectively; 

P:  tunneling  probability. 

Note  that  (2.1)  simplifies  to  the  commonly  used  direct  tunneling  current 
expression  [MajOO,  Sha99,  She72,  Dep95]  when  mta*  and  mtb*  are  equal, 
e.g.,  in  the  case  of  conduction-band  electron  tunneling. 

The  model  in  (2.1)  assumes  conservation  of  Et  as  well  the  total 
energy  (E  = Et  + Ex)  in  the  (elastic)  tunneling  process,  which,  based  on  the 
WKB  approximation  for  the  matrix  element  for  the  transition  [Har61]  and 
a one-band  parabolic  dispersion  relation  involving  an  effective  mass 
(mox*)  in  the  Si02  [Dep95],  yields 
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(2.2) 


where,  from  Figure  2.3,  EC(ox)(x)  = Evs  + (|>bv  - q£0Xx  is  the  conduction 
band-edge  energy  in  the  oxide;  <|)bv  = 4.2eV  and  £ox  = V|/ox/tox  is  the 
(constant)  electric  field  in  the  oxide.  However,  (2.1)  and  (2.2)  do  not 
explicitly  account  for  the  conservation  of  transverse  momentum  (via 
phonons)  in  this  case  of  an  indirect  valence  band-conduction  band 
transition  [Mol64].  This  accounting  decreases  P dramatically,  but  the 
functional  dependence  on  £ox  is  almost  identical  [Mol64]  to  (2.2),  but  with 
mox*  being  a larger  effective  mass  [Sze81].  Several  papers,  e.g.,  [Dep95], 
report  different  electron  effective  masses  in  the  oxide,  corresponding  to 
conduction  band-conduction  band  tunneling  and  depending  on  tox  and 
model  assumptions.  We  infer  from  measurements  of  the  valence-band 
tunneling  current  (discussed  later)  mox*  = 0.36m0,  where  m0  is  the  free 
electron  mass. 


verified  numerically)  that  we  can  assume  (with  <1%  error)  fb  = 0.  Thus,  for 
the  tunneling  component  of  interest,  (2.1)  simplifies  to 


where  mv*  (=  mta*)  is  the  transverse  DOS  effective  mass  for  electrons  in 
the  silicon  valence  band.  We  assume  mv*=0.65m0,  which  follows  from 
summing  the  heavy  (mvb*=0.49m0)  and  light  (mvi*=0.16m0)  hole  masses 


Referring  to  Figure  3(c),  we  see  that  fa  = 1,  and  (as  we  have 


(2.3) 
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[Sze81]  in  the  valence  subbands  to  properly  account  for  the  composite  2D 
DOS.  Note  in  Figure  3(c)  that  the  total  energy  (E)  of  a tunneling  electron 
contributing  to  the  valence-band  component  of  interest  here  must  be 
between  ECG  and  Evs;  the  tunneling  to  states  in  the  gate  below  ECG  is 
negligible.  We  can  thus  define  the  Ex  limits  of  integration  in  (2.3)  as  Evs 
- Et  and  ECG  - Et;  and  Et(max)  is  Evs  - Ex.  However,  Ex  and  Et  are  correlated 
in  (2.3),  and  hence  we  need  a simplifying  assumption  to  do  the  integration. 
We  assume  that  in  the  valence  band,  electron  kinetic  energy  is  negligible 
compared  to  the  potential  energy,  or  Et  = 0.  Therefore,  (2.3)  is 
approximated  as 

jvbe  - ^cjrE'pdE>dE>=  <2-4> 


Ecg*  0 


K 


JCG 


To  evaluate  (2.4),  (2.2)  is  first  integrated  for  specific  Ex: 


P = exp 


AtoxJ2mox*^bv  +EVS~Ex)  EVS~Ex-(lVox) 

3 qh  V|/ox 


(2.5) 


With  (2.5)  substituted  in  (2.4),  we  note  that  (Evs  - Ex)  is  a much  weaker 
function  than  the  exponential  function.  Therefore,  it  can  be  removed  from 
the  integration,  provided  a correction  factor  is  inserted.  Numerical 
solution  show  that  a factor  of  1/2  is  acceptable;  thus, 


J 


VBE- 


4nmv*q(Evs-ECG ) 
h3 
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4toxJ2mox*(<lhv  + EVS 

3qh 


Ex)  -(Q$bv  + Eys~  Ex-g\\fox) 

Vox 


3-, 

2 


dEx, 


(2.6) 


which  we  rewrite  as 


dvBE=2^Evs~ECG) 


r^VS 

L 1 

J EcG 


exp 


- Btr 


3 3. 

(q§bv  + Evs  - Ex )2  - (qtybv  + Evs  - Ex  -qyox) 

c Vox 


dEv  (2.7) 


Anm*q 

A = 5 — , n = 


HT  3(lh 

Integrating  (2.7),  and  noting  that  for  inversion  conditions 


Btr 


Vo 


Qhv  -(qhv-Wox) 


» i 


Btr 


Vc 


(<7<l>bi>  + Evs-Ecg)  ~(q<\>bv  + EyS-  ECG-q\\iox) 


»1, 


we  obtain 


J 


VBE 


A ^ ox  p 

3BV„  ‘ 


F1exp 


- Btr 


Q<\>bv  -(g^bo-gVox) 

c \l/ 

T OX 
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tH 

q<\>bv2-(q<bbv~qVox)2 

,F2  = 

q<\>bv'2-(q<bbv'-qVox)2 

Ed  = Eys~E  CG  ’ Q§bv  - Q§bv  + Ed- 

From  Figure  2.3,  for  EF  = Ec  in  the  highly  doped  poly-gate  region,  we  can 
write 


Ed  ~ Evs-Ecg-  q 


E 

VoB-Vsb-Vp-^j;  + <t>^]; 


(2.9) 


Eg  is  the  silicon  bandgap  and  <])F  is  the  Fermi  potential  in  the  substrate. 
We  note  that  JVBE  in  (2.8)  must  be  zero  when  Ed  is  zero  or  negative.  To 
validate  the  tunneling  current  model  then,  Ed  in  (2.8)  is  replaced  by  the 
smoothing  function, 


— - ln{  1 + exp[SyBE(.Ed/<7)]} 

Ed/q=  


S 


(2.10) 


VBE 


where  SVBE  is  a constant.  For  Ed  > 0,  Ed  = Ed  \ but  when  Ed  decreases  to 
zero  and  becomes  negative,  Ed  approaches  zero,  which  properly  forces 
JVBe  to  zero. 

For  a preliminary  check  of  the  JyBE  model  defined  by  (2.8M2.10), 
we  combine  it  with  a numerical  solution  of  the  basic  MOS  equation 
[Tau98], 

= VS6  + Vo*  + > (2-11} 

coupled  to  the  1-D  (in  x)  Poisson  equation;  OpS  is  the  poly-Si/Si  work- 
function  difference,  which  is  approximated  by  assuming  that  the  Fermi 
level  in  the  bulk  gate  region  is  pinned  at  EC(gate)  due  to  the  high  doping 
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density,  and  yp  is  the  potential  drop  in  the  poly-Si  gate,  which  is 
characterized  based  on  the  depletion  approximation  [Tau98].  We  show  in 
Figure  2.4  model-predicted  JVBE  f°r  three  different  values  of  tox  (1.65nm, 
1.79nm,  and  2.18nm,  which  were  determined  via  ellipsometer 
measurements),  with  measured  gate-to-body  tunneling  currents 
superimposed.  The  measured  data  were  obtained  from  bulk-silicon 
nMOSFETs  via  a carrier-separation  experiment  as  described  in  [Shi98]. 
The  gate  was  positively  biased  (VGB),  while  the  source/drain  and  the 
substrate  were  grounded.  The  source/drain  current  (Js/d)  thereby 
measured  is  component  ® in  Fig.  1,  while  the  measured  substrate  current 
(JB)  is  component  ®.  The  model  predictions,  with  mox*  tuned  to  0.38m0 
and  SVBE=14.0,  are  in  g°od  agreement  with  the  measured  data  for 
relatively  low  VGB,  but  underestimate  the  gate-to-body  tunneling  current 
for  high  VGB,  especially  for  thicker  gate  oxide. 

2.3.2  Scaling  Effects 

The  measured  data  in  Figure  2.4  were  taken  from  scaled 
nMOSFETs.  The  discrepancies  between  the  data  and  the  JyBE  model 
predictions  thus  could  possibly  be  due  to  scaled-device  effects  that  we 
heretofore  ignored.  In  this  section,  we  consider  three  such  effects:  electron 
degeneracy  in  the  gate  and  electron  energy  quantization  and  exchange 
energy  in  the  inversion  layer.  The  noted  degeneracy,  due  to  the  high  gate 
doping  density,  results  in  Ep  being  above  EG  in  the  gate,  with  the 
difference  q = (EF  - Ec)/(kBT)  defined  by  Fermi-Dirac  statistics  [Sze81]: 


Gate-to-Body  Current  Density  (A/cm2) 
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Figure  2.4  Measured  (points)  and  modeled  (curves)  gate-to-body 
tunneling  current-density  (Jvbe)  characteristics  in 
nMOSFETs  with  varying  gate-oxide  thickness.  The  gate  area 
is  lOxlOmm2,  the  substrate  doping  density  is  1.7xl017cm'3, 
and  the  n+-poly-Si  gate  doping  density  is  approximately 
lxl020cm"3;  VBg  = V^g  = 0,  mox*  = 0.38m0. 
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Np  = NcF1/2(T\) 


(2.12) 


where  Np  is  the  n+poly-gate  doping  density  and  Nc  is  the  effective  density 
of  states  in  the  conduction  band.  (For  the  p+  poly-Si  gate  on  the  pMOSFET, 
Nc  is  replaced  by  Nv,  the  effective  density  of  states  in  the  valence  band.) 
We  use  Nilsson’s  approximation  [Nil78]  for  r\: 


In  this  approximation,  Ep  - Eq  is  57.0meV  when  Np>  is  1020cm  3,  which  is 
significant.  Thus,  the  substrate-gate  overlapped  energy  for  tunneling, 
(2.9),  should  be  replaced  by 


where  Vx  = kBT/q  is  the  thermal  voltage. 

The  noted  quantization  effect,  for  a specified  VqB,  increases  \|/sb 
[Tau98],  and  hence  can  affect  JyBE  [Yan99].  In  the  UFPDB  model  [Pel02], 
the  quantization  effect  is  implemented  [ChiOl]  employing  van  Dort’s 
model  [Van94],  where  the  effect  is  modeled  by  calculating  an  effective 
intrinsic  carrier  density  corresponding  to  the  bandgap  widening  AEg(QM) 
implied  by  the  creation  of  2-D  electron  subbands  above  Ec.  However,  the 
noted  exchange  energy  lowers  the  energy  of  the  subbands,  thereby 
effectively  narrowing  the  bandgap.  We  use  Stern’s  model  [Ste73]  to 
account  for  this  narrowing: 


(2.13) 


Ed~  Evs~Ecg  ~ Q Vgb~ 


(2.14) 
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l 


AEg(EX ) - 


(2.15) 


where  Nj  is  the  areal  electron  density  (Nj  = -Q/q  where  Qj  is  the  inversion 
charge  density  approximated  as  -£giExf0  [Tau98];  Exf0  is  the  surface 
electric  field),  nv  is  the  multiplicity  of  the  valleys  from  which  the  lowest 
subband  arises  (nv  = 2 for  inversion  electrons,  nv  = 1 for  inversion  holes), 
and  k is  the  effective  static  dielectric  constant,  which  is  (kgj  + kox)/2  in  the 
extreme  2D  electron-limit  case.  For  strong  inversion,  we  thus  modify  the 
bandgap  widening  to  account  for  the  exchange-energy  narrowing: 


Now,  for  specified  Vqq  (=  Vqb  with  Vgg  = 0)  with  (2.16),  the 
model-defined  ysb  and  Qj  (from  UFPDB  discussed  in  Sec.  2.5)  are  used 
with  (2.8)  and  (2.14)  to  define  JyBE-  For  the  tox  = 1.65nm  nMOSFET,  we 
show  in  Figure  2.5  the  model-predicted  Jvbe  with  and  without  the  three 
scaled-device  effects,  compared  with  the  measured  substrate  current.  Note 
that  when  the  quantization  effect  is  accounted  for,  both  the  degeneracy 
and  exchange-energy  effects  must  be  modeled  in  order  to  get  reasonable 
agreement  with  the  data.  This  agreement  is  shown  in  Figure  2.6  for  all 
values  of  tox.  The  better  agreement  between  the  measured  and  model- 
predicted  Jvbe  shows  the  importance  of  the  noted  scaled-device  effects. 
For  all  three  devices,  a tuned  value  of  mox*  = 0.36m0,  which  is  lower  than 


A Eg  = A E 


' g{QM)~  AE g{EX)  • 


(2.16) 


the  value  inferred  previously  without  the  scaled-device  effects,  was  used. 
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VGB(V) 


Figure  2.5  Modeled  valence-band  electron  tunneling  current  density  of 
the  tox=1.65nm  nMOSFET,  showing  the  effects  of 
quantization  (QM),  gate  degeneracy  (DEG),  and  channel 
exchange  energy  (EX);  Vgg  = Vds  = 0 and  mox*  = 0.36m0.  The 
measured  gate-to-body  current  is  also  shown. 
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Figure  2.6  Measured  and  modeled  gate-to-body  tunneling  current- 
density  (Jvbe  with  scaled-device  effects)  characteristics  of 
nMOSFETs  with  varying  gate-oxide  thickness;  VBs  = Vos  = 
0,  mox*  = 0.36m0. 
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The  agreement  is  good  for  low  Vqb>  but  large  discrepancies  remain  for 
high  VGB. 

2.3.3  Impact-Ionization  Effects 

The  remaining  model-measurement  discrepancies  evident  in 
Figure  2.6  are  similar  to  those  reported  in  [Sha99]  for  thin  tox,  although 
the  model-predicted  curves  in  [Sha99]  were  tuned  for  better  fits  to  the 
data  at  high  VGB.  The  discrepancies  in  [Sha99]  were  attributed  to 
limitations  of  the  WKB  approximation  for  ultrathin  oxides.  However,  we 
believe  electron-hole  pair  generation  by  impact  ionization  in  the  gate  and 
subsequent  valence-band  hole  tunneling  underlie  the  large  discrepancies 
in  Figure  2.6  and  define  the  dominant  gate-to-body  direct  current 
component  at  high  VqB. 

As  illustrated  in  Figure  2.1(a),  the  conduction-band  electrons 
tunneling  from  the  substrate,  for  high  VGB,  become  hot  and  have  sufficient 
kinetic  energy  to  drive  significant  impact  ionization  in  the  gate.  A portion 
of  the  hot  electrons  will  lose  their  energy  via  electron-electron 
interactions,  but  we  presume  some  surface  depletion  in  the  poly-Si  gate 
that  reduces  the  electron  density  enough  to  enable  a significant  fraction  of 
the  electrons  to  lose  their  energy  via  impact  ionization.  Holes  in  the 
valence  band  generated  by  the  impact  ionization  can  then  tunnel  to  the 
substrate  valence  band,  thus  significantly  augmenting  component  © in 
Figure  2.1(a).  To  check  this  theory,  we  assume  a steady-state  condition 
where  a hole  density  is  established  in  the  gate,  and  the  hole  tunneling 
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current  is  defined  directly  by  the  rate  of  hole  generation  via  the  impact 
ionization.  This  assumption  presumes  that  the  tunneling  process  is  much 
faster  than  the  hole-electron  recombination  process  in  the  gate,  as  well  as 
that  the  hole  tunneling  without  the  impact  ionization  is  negligible.  Thus, 
with  reference  to  Figure  2.3(c),  we  express  the  gate-substrate  hole 
tunneling  current  as 


J VBH  - C {M-l)dJCBE  (2.17) 

J &CS 

where  J’cbe  is  the  mentioned  fraction  of  the  conduction  band  (channel)- 
gate  electron  tunneling  current  Jcbe  (component  ® in  Figure  2.1(a)),  Eqq 
is  the  conduction-band  edge  energy  at  the  silicon  surface,  and  (M-l)  is  the 
impact-ionization  multiplication  factor,  which  depends  on  the  kinetic 
energy  (E  - E^q)  of  the  electrons  in  the  gate;  E is  the  total  energy  of  the 
tunneling  electrons.  We  express  the  kinetic  energy  in  terms  of  the  electron 
temperature  [Sah91]: 

E-Ec  G = | kBTe.  (2.18) 


For  a first-order  characterization  of  the  impact  ionization,  we  write 


(Af-  l)  = aXe  (2.19) 

where  a,  dependent  on  Te,  is  the  ionization  rate  and  A.e  is  the  energy 
relaxation  length.  For  quasi-steady-state  conditions  [Slo91], 


oc(Te)  = a0  exp 


rjoi 
. y(Te-T)_ 


(2.20) 
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where  y = 5kB/(2qA,e),  a0  and  p0  are  constants,  and  T is  the  lattice 
temperature. 

With  (2.18M2.20),  (2.17)  with  models  for  Jcbe  anc*  J’cBE  can  be 
integrated  over  the  conduction  band  to  characterize  JyBH-  To  corroborate 
our  theory,  however,  we  assume  J’cBE  ~ ^CBE  and  use  measurements  of 
the  conduction-band  electron  tunneling,  which  we  assume  occurs 
predominantly  at  or  near  Eq§  + AEg  where  AEg  is  a effective  bandgap 
widening  in  (2.16).  Then,  the  integral  in  (2.17)  is  expressed  as 


r (M -l)dJCBE  = (M -l)JCBE(me 

J^CS 


as) 


where  (M-l)  is  evaluated  at  E = E^s  + AEg,  yielding 


(2.21) 


Jvbh  - aKJcBE(meas)  (2.22) 

with  a evaluated  from  (2.20)  at  Te  = 2(Eqs  + AEg  - EcG)/(3kB),  which 
follows  from  (2.18);  the  presumed  kinetic  energy  can  be  expressed  in  terms 
of  the  gate  voltage  as  evident  in  Figure  2.3(c): 


Ecs  + A Eg  - Ecg=  Ed  + Eg  + A Eg  (2.23) 

where  Ed  = EVs  - ECG  from  (2.14).  In  (2.23),  we  have  implicitly  assumed 
that  the  transverse  kinetic  energy  is  negligible  since  it  is  not  included  in 
AEg.  (For  high  Vds,  the  inversion  electrons  can  have  high  transverse 
kinetic  energy  near  the  drain  due  to  high  electric  field  in  channel. 
However,  the  gate-to-body  direct  tunneling  current  is  negligible  in  this 
region,  because  of  relatively  low  Ed  and  eox,  and  hence  the  noted 
assumption  is  innocuous.)  Figure  2.7  shows  the  kinetic  energy  in  (2.23) 
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Figure  2.7  Simulated  kinetic  energy  of  the  tunneling  electron  from  the 
conduction  band  of  the  substrate.  The  onset  voltage  for 
impact-ionization  is  indicated  for  varying  oxide  thickness. 
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versus  VGB,  derived  via  UFPDB,  for  the  three  values  of  tox  in  Figure  2.6. 
Our  first-order  characterization  of  Jvbh  in  (2.22)  is  applicable  when  the 
kinetic  energy  is  greater  than  the  physical  onset  value  for  impact 
ionization.  We  define  this  onset  via  a classical  ballistic  model  governed  by 
conservation  of  both  energy  and  momentum  [Sah91];  it  estimates  the 
onset  energy  as  3Eg/2.  Figure  2.7  shows  that  the  onset  VGB  increases  with 
decreasing  tox  as  implied  by  Figure  2.6. 

For  the  tox  = 1.65nm  device  of  Figure  2.6,  we  plot  in  Figure  2.8 
JVBh(VGb)  evaluated  from  (2.22)  for  VGB  greater  than  the  onset  voltage 
(2.2V,  as  inferred  from  Figure  2.7).  For  the  evaluations,  we  used  the 
measured  JcBE(meas)(^GB)  (=  jS/d)>  a^so  plotted  in  Figure  2.8,  and  we 
assumed  a0  = 2.45xl06cm'1  and  Xe  = 65nm  as  given  in  [Slo91],  and  we 
tuned  po  = 2.98xl06V-cm'1,  which  is  larger  than,  but  comparable  with  the 
value  given  in  [Slo91].  The  larger  p0  could  reflect  some  energy  dissipation 
via  electron-electron  interaction  in  the  gate.  For  comparison,  we  include 
in  Figure  2.8  the  measured  gate-to-body  direct  tunneling  current  (JB)  and 
JVBe(VGb)  from  (2.8),  which  are  in  Figure  2.6,  and  the  predicted  sum 

(JvB(total))  of  JVBE  and  JVBH-  Note  how  wel1  JVB(total)(VGB)  matches  JB  for 
all  VGb  in  the  inversion  region,  and  that  JyBH>  driven  by  the  impact 
ionization,  is  the  dominant  component  for  high  VGb-  Note  also  the  steeper 

slope  of  JVBh(VGb)>  relative  to  that  of  JVBe(VGb)  and  JCBE(meas)(VGB)-  The 
VGB  dependence  of  (M-l)  in  (2.21)  causes  the  steeper  slope,  which  cannot 
be  explained  by  possible  uncertainties  in  tox,  mv*,  or  mox*  in  the 
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Figure  2.8  Measured  and  modeled  gate-to-body  current-density 
characteristics  of  an  nMOSFET,  including  the  two  predicted 
valence-band  components;  Vgg  = 0.  The  impact-ionization 
effect,  driven  by  the  conduction-band  electron  tunneling 
component,  the  measured  characteristic  of  which  is  included, 
is  illustrated. 
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Jvbe(VGb)  model  while  maintaining  good  agreement  with  the  measured 
JB  for  lower  VGB. 

Using  the  same  impact-ionization  model  with  the  same 
parameter  values,  we  evaluated  JygH  and  JvB(total)  f°r  other  devices 
in  Figure  2.6  having  different  tox;  the  model  predictions  and  the  measured 
JB  are  plotted  versus  VGB  in  Figure  2.9.  We  used  the  onset  voltages  for 
impact  ionization  evident  in  Figure  2.7.  The  good  agreement  between  the 
measured  and  model-predicted  Jg  at  high  VGB  in  Figure  2.9  supports  our 
theory  of  a significant  impact-ionization  effect.  (Note  that  a smoothing 
function  similar  to  that  in  (2.10)  would  have  to  be  used  in  the  model  to 
effect  the  proper  turn-on  of  the  impact  ionization.) 


We  next  consider  the  gate-to-body  tunneling  for  VGB  < 0 as 
reflected  by  the  energy-band  diagrams  in  Figure  2.10.  For  accumulation, 
the  general  tunneling  current  expression  (2.1),  with  mta*  = mtb*  = mc*  in 
the  conduction  band,  describes  the  predominant  conduction-band  electron 
tunneling,  i.e.,  component  ® in  Figure  2.1(b),  as 


where  Et(max)  can  be  virtually  infinite;  me*  is  the  transverse  DOS  effective 
mass  for  electrons  in  the  (100)  silicon  conduction  band  (mc*  = 0.19m0). 
Noting  in  Figure  2.10(b)  that  the  minimum  Ex  must  be  ECG,  and 


2.4  Depletion/Accumulation-Region  Model 


(2.24) 
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Figure  2.9  Measured  and  modeled  gate-to-body  current-density 
characteristics  of  nMOSFETs  with  varying  gate-oxide 
thickness;  VBS  = VDS  = 0>  mox=0.36m0,  and  SVBE=13.5. 
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(a)  Depletion 


A 


(b)  Accumulation 


Figure  2.10  Energy  bands  showing  the  conduction-band  electron 
tunneling  for  the  depletion  and  accumulation  cases  in  the 
nMOSFET.  The  valence-band  hole  tunneling  is  also  shown  as 
a dominant  component  for  strong  accumulation  in  (b). 
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integrating  (2.24),  with  the  Fermi-Dirac  distribution  functions  f(E  = Et  + 
Ex),  over  Et,  we  get 


r Er  — E 

1 + 


CBE 


47 zqmc 


kTf  P(EX)  In 


dE . 


(2.25) 


A classical  WKB  approximation  and  a one-band  parabolic  dispersion 
relation  was  used  to  calculate  the  tunneling  probability: 


P = 


(2.26) 


In  Figure  2.10, 

Ecox(x)=  §bc  + ecg  + QEoxx  (2.27) 

where  eox  is  the  electrical  field  strength  in  the  oxide  layer  and  §bc  is  the 
tunnel  barrier  height.  Thus, 

k = — 72mox*(<t)fec  + ecg  ~Ex  + Q£oxx)  (2.28) 

and 


P = exp 


4£, 


3 qH 


+ ECG 


3 3- 

Ex  + g^Qx)2  ~ + ECG  ~ Ex ) 

Vox 


(2.29) 


where  \j/ox  = eoxtox . 

In  Figure  2.10(b), 

Ex~EFG~  (Ex  ~ ECg)  ~ (2.30) 


Ex~efs  ~ (Ex~efg)-QVgb 


(2.31) 
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where  r\  (=  (EpG  - EGG)  / (kgT))  is  the  degeneracy  due  to  the  high  gate 
doping  density  in  (2.13).  The  integration  in  (2.26)  is  done  numerically 
(trapezoidal)  with  a finite  upper-Ex  limit  (EGG  + 0.3eV),  partitioning  the 
integration  range  into  50  equal  segments. 

As  illustrated  in  Figure  2.10(a),  a finite  Jcbe  actually  flows  prior 
to  the  onset  of  strong  accumulation.  The  tunneling  current  flows  to  the 
channel,  body,  or  both,  depending  on  VGg.  Therefore,  for  modeling  the 
gate-to-body  component,  we  need  the  onset  condition  of  gate-to-body 
current  and  it  is  defined  by  following  smoothing  function,  which  properly 
forces  the  gate-to-body  current  to  zero  when  the  surface  potenital  is  above 
the  some  fraction  of  <t>p  and  a significant  channel  exists: 


J CBE  - JcBEX 


1 + exp 


\|/s6-Fx<t>F 


}CBE 


(2.32) 


With  (2.32),  the  formalism  of  the  tunneling  current  in  the  accumulation 
bias  region,  (2.25),  is  applicable  to  calculate  the  tunneling  current  in  the 
depletion  bias  region  when  the  lower  limit  of  the  integral  is  EGg  instead  of 
Ecg.  However,  we  need  to  smooth  the  change  of  the  lower-energy  limit  in 
each  bias  region.  The  following  smoothing  function  is  used: 


ecg=  Ecs  + Q lni  1 + exP 


' FB 


ECG~ECS 


FB 


(2.33) 


where  Spg  is  a constant.  For  EGG  < EGg,  Eqq  = EGg;  but  when  EqG 
increases  above  Ecs,  ECG  = EGG,  which  properly  change  the  lower  limit. 
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For  a preliminary  check  of  our  model  for  Jcbe  in  (2.25),  we 
combine  it  with  a numerical  solution  of  the  basic  MOS  equation  (2.11), 
without  i|/p  (i.e.,  without  poly  depletion  for  accumulation)  but  with  the 
degeneracy  effect  in  (2.13)  included  in  <hPS,  coupled  to  the  1-D  (in  x) 
Poisson  equation.  We  show  in  Figure  2.11  the  model-predicted  Jcbe  for 
three  different  values  of  tox  (1.65nm,  1.79nm,  and  2.18nm),  with  measured 
gate-to-body  tunneling  currents  superimposed.  The  model  predictions, 
with  mox*=0.36m0,  the  same  value  finally  used  for  the  JVBE  model,  are  in 
good  agreement  with  the  measured  data  for  relatively  low  | VGB  | in  the 
depletion/accumulation  region,  but  underestimate  the  gate-to-body 
tunneling  current  (magnitude)  for  high  | VGB  | . This  difference  is  likely 
due  to  the  valence-band  hole  tunneling,  depicted  in  Figure  2.10(b),  which 
becomes  significant  for  strong  accumulation  of  holes.  Since  normal 
operation  of  scaled  CMOS  devices  would  avoid  this  region,  we  do  not 
include  this  component  of  tunneling  current  in  our  model.  For  the  tox  = 
2.18nm  device,  the  error  seems  to  be  prevalent  even  for  lower  | VGB  | , but 
we  note  in  Figure  2.11  that  this  apparent  discrepancy  is  probably  due  to 
the  low  values  of  JGBe  that  are  difficult  to  separate  out  from  the  noise  in 
the  measurement  system. 

2.5  Implementation  of  Gate-to-Body  Direct  Tunneling  Current  in  UFPDB 
The  complete  gate-to-body  direct  current  model  was 
implemented  in  UFPDB  [Pel02,  ChiOl]  for  reliable  prediction  of  floating- 
body  effects  in  scaled  PD/SOI  CMOS  devices  and  circuits.  We  overview  in 
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Figure  2.11  Measured  (points)  and  modeled  (curves)  gate-to-body 
current-density  (magnitude  of  Jcbe  < 0)  characteristics  of 
nMOSFETs  with  varying  gate-oxide  thickness;  VDS  = 0,  mox* 
= 0.36m0,  Sqbe=0  04,  F=0.5,  and  Spg=10.0. 
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this  section  the  implementations  of  Ivbe^GS»  Vds)  and  IcBE^GS>  VBS, 
VDs)- 

In  Sec.  2.3,  Ivbe(VGb)  = ^VBE^GS  = VGB>  VBS  = °>  VDS  = °)  was 
developed.  For  arbitrary  bias,  we  write  VGB  = VGS  - VBS  and  \j/sb  = \|/sf-  VBS, 
where  \|/sf,  used  in  UFPDB,  is  the  absolute  surface  potential  in  the 
channel,  referenced  to  that  in  the  hypothetical  unbiased  body  region.  We 
note  from  (2.14)  that  we  are  hence  using  the  Ivbe(Vgb)  formalism  directly 
with  VGB  and  \|/sb  replaced  by  VGS  and  \j/sf,  respectively. 

We  extend  the  model  to  include  VBg  dependence  as  follows.  The 
total  valence-band  electron  tunneling  current  is  defined  by  integrating 
JVBe  along  the  channel: 

^VBE  ~ ^VBE^sf^y  (2.34) 


where  Leff  is  the  effective  channel  length  and  W is  the  channel  width.  For 
strong  inversion, 

Vsf(y)  = vSfo  + v(y)  (2.35) 

where  \|/sf0  = \j/sf(y  = 0)  is  the  surface  potential  at  the  source  end  of  the 
channel  and  V(y)  is  the  quasi-Fermi  potential  for  the  inversion  electrons. 
Assuming  that  the  electron  channel  current  is  constant  in  y,  we 
approximate  V(y)  as  [Tau98] 


V(y)  = 


Vors-Vt 

m 


Vofs-Vt 

m 


2 y ( VGfS  ~ Vt 


JefP 


m 


V°S+  t 


-V 


l 

-i2 


eff 


DS 


(2.36) 
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where  m = 1 + Cd/Cox;  Cox  is  the  gate-oxide  capacitance,  Cd  is  the  body 
depletion  capacitance,  and  Vt  is  a conventionally  defined  threshold 
voltage.  Since  UFPDB  does  not  involve  Vt  directly,  we  assume  in  (2.36) 
that 


^GfS~  V* 
m 


= V 


DS(sat) 


(2.37) 


where  VDS(sat)  is  the  saturation  drain  voltage  defined  in  UFPDB  [ChiOl] . 
Further,  to  smoothly  account  for  saturation,  we  replace  Vpg  in  (2.36)  by 
VDSx,  a UFPDB-defined  smoothing  function  that  equals  VDS  for  low  VDS, 
but  is  limited  to  VDg(sat)  at  high  VDg. 

To  complete  the  characterization  of  in  Ivbe  in  (2.34),  we  need  to 
characterize  yp(y)  in  (2.14).  This  is  done  in  UFPDB  [ChiOl]  via  the 
depletion  approximation  in  the  poly-Si  gate,  which  yields  the  depletion 
charge  in  terms  of  \j/p: 

Qp  = qNpxp  = j2eSiqNp\\ip  . (2.38) 

Applying  Gauss’  law  to  the  gate-oxide  interface  gives 


- ~t  “ - Cox(VGS-®PS-Vsf-  Vp) 


(2.39) 


which,  with  (2.38),  yields 


VpCy)  - -~2'((lNPS'Si  + ^ox(VGS~^>PS~x\'sf') 
Cox 


-jQNP£si[QNpESi  + 2C0X(Vg/-s  - - \|/s^)])  ,(2.40) 


with  ysf{y)  given  by  (2.35). 
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Now,  for  a specified  bias  (VGs»  VDS),  the  integrand  JVBE  in  (2.34) 
is  evaluated  via  (2.8),  (2.14),  and  (2.35)-(2.40),  using  the  strong-inversion 
VsfO  determined  in  UFPDB  [ChiOl].  To  properly  incorporate  the 
quantization  effect  and  the  exchange  energy  into  the  vj/sf0  calculation,  a 
few  Newton-Raphson  iterations  are  required  since  those  effects  and  ysf-0 
are  correlated.  For  the  first  iteration,  \|/sf0  is  a strong-inversion  guess 
based  only  on  the  device  structure.  The  quantization  effect,  dependent  on 
AEg(QM)  [Van94],  and  the  exchange  energy  dependent  on  Qj  in  (2.15)  are 
calculated  using  V|/sfQ.  For  the  subsequent  iterations,  an  effective  intrinsic 
carrier  density  corresponding  to  the  bandgap  widening  AEg  in  (2.16)  yields 
the  new  \(/sfQ  iterate.  Generally,  four  iterations  are  used. 

The  integration  in  (2.34)  is  then  done  numerically  (trapezoidal 
method),  partitioning  Leff  into  10  equal  segments.  Note  that  the 
temperature  dependence  of  IVBE  is  defined  implicitly  via  ysf(T),  and 
explicitly  via  NC(T),  Eg(T),  fF(T),  and  VT  in  (14).  Figure  2.12  shows  the 
simulated  IyBE  characteristics  of  an  nMOSFET  with  varying  VGS  and  VDS. 
Note  how  the  integrated  tunneling  current  decreases  substantively  with 
increasing  VGg  due  mainly  to  (Eyg  - EGG)  decreasing  along  the  channel. 

In  the  accumulation  region,  the  integrated  IGbe  can  be  expressed 
simply  as 

^CBE  - WLeffJ GBE^sb*  (2.41) 

since  the  surface  potential  is  constant  along  the  channel.  In  (2.41),  q/sb, 
which  is  calculated  by  subtracting  VBg  from  the  UFPDB-defined  [ChiOl] 
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Figure  2.12  UFPDB-simulated  valence-band  electron  tunneling  current 
(IVBe)  characteristics  of  an  nMOSFET  for  varying  Vqs  and 
VDs;  Leff  = W = 10mm,  tox  = 1.65nm,  mox*  = 0.36m0. 
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Vsf<y  = 0).  is  smoothed  to  zero  when  it  tends  to  be  less  than  zero  in  the 
accumulation  region.  The  total  gate-to-body  current,  IGB,  in  UFPDB  is 
defined  by  summing  IVBE  > 0 in  (2.34)  and  ICBE  < 0 in  (2.41):  IGB  = IVBE  + 
icBE- 

Since  the  IGB(VGg,  VBS,  VGg)  formalism  is  complex,  we  use 
difference  approximations,  based  on  AVGS  and  AVDS  perturbations,  in 
UFPDB  to  calculate  the  derivatives  needed  for  the  SPICE  (Newton- 
Raphson)  nodal  analysis.  The  perturbed  values  of  IyBE  follow  directly  from 
(2.34),  with  \|/sf0  virtually  fixed  because  of  the  strong-inversion  condition 
(VGg  > V^g  [ChiOl]).  The  electron  effective  mass  in  the  oxide  layer  (mox*) 
and  the  smoothing-function  constant  (SVBE,  SCBE,  and  F)  are  the  new 
model  parameters  needed.  Further,  mox*  is  used  as  a flag  to  turn-on  the 
IGB  model:  if  mox*=0,  then  the  IGB  formalism  is  bypassed  (although  the 
accounting  for  the  exchange  energy  via  (2.16)  is  permanent  in  UFPDB-2.0 
for  VGfg  > VTS);  SyBE  and  SCBE  are  used  as  flags  to  turn-on  the  IVBE  and 
IGBE  models,  respectively. 

For  the  pMOSFETs,  with  same  values  of  tox  as  the  nMOSFETs 
previously  modeled,  we  show  in  Figure  2.13  the  UFPDB-predicted  JyBE- 
VGB  characteristics.  The  model  parameter  mox  is  tuned  as  0.37m0,  versus 
0.36m0  for  the  nMOSFETs.  Since  the  tunneling  processes  in  both  devices 
are  identical,  the  slightly  different  values  of  mox  possibly  reflects  slight 
differences  in  the  device  structure,  e.g.,  in  tox.  The  model-measurement 
agreement  in  Figure  2.13  is  very  good  for  low  and  moderate  VGB,  but  large 
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Figure  2.13  Measured  and  modeled  (Jvbe  with  scaled-device  effects) 
gate-to-body  current  characteristics  of  pMOSFETs;  Vgg  = 
VDS  = 0,  mox  = 0.37m0. 
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discrepancies  remain  for  high  VGB  as  we  found  for  the  nMOSFETs.  Thus, 
the  impact-ionization  effect  at  high  VGB  appears  to  occur  in  the 
pMOSFETs,  as  depicted  in  Figure  2.2(a),  as  well  as  in  the  nMOSFETs.  As 
shown  in  Figure  2.13,  the  onset  voltages  of  impact-ionization  is  slightly 
lower  than  those  of  nMOSFETs  because  the  kinetic  energy  of  tunneling 
electron  in  pMOSFETs  expressed  by  (2.23)  is  bigger  than  that  in 
nMOSFETs  due  to  the  higher  degeneracy  in  p+poly  gate. 

In  Figure  2.14,  we  show  the  the  UFPDB-predicted  JyBH'^GB 
characteristics  for  the  pMOSFETs.  We  don’t  include  the  conduction-band 
electron  tunneling,  which  is  dominant  in  deep  accumulation  region,  since 
typical  operation  voltage  in  this  scaled  technology  does  not  go  to  the 
region.  The  valence-band  hole  tunneling  formalism  is  identical  to  the 
conduction-band  electron  tunneling  in  (2.25)  except  the  transverse  DOS 
effective  mass  for  electrons  in  the  silicon  conduction-band  (mc*=0.19m0), 
which  should  be  replaced  by  the  transverse  DOS  effective  mass  for 
electrons  in  the  silicon  valence  band  (mv*=0.65mo).  The  model  predictions, 
with  mox*=0.37mo,  which  is  the  same  value  as  in  JyBE  modeling, 
£>CBE=0-045,  f=o.8,  and  Sjrg=15.0,  are  in  good  agreement  with  the 
measured  data  for  depletion/accumualtion-region  with  relatively  low  Vqb- 

Continuous  scaling  of  the  Si02  gate-dielectric  thickness  (tox)  has 
resulted  in  significant  direct  gate  tunneling  current,  as  noted  herein,  not 
to  mention  the  strong  tendency  for  penetration  of  the  poly-gate  dopant 
into  the  silicon-channel  region.  Alternative  gate  dielectrics,  having  higher 
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Figure  2.14  Measured  and  modeled  (Jvbh)  gate-to-body  current-density 
characteristics  of  pMOSFETs  with  varying  gate-oxide 
thickness;  Vbs  = VDg  = 0,  mox  = 0.37m0,  S^be  = 0.045,  F = 0.5, 
and  SpB  = 15.0. 
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dielectric  constants  and  hence  enabling  thicker  tox,  yet  high  gate 
capacitance  and  drive  current,  are  thus  being  examined  for  ultra-scaled 
CMOS  [Shi99] . Therefore,  we  further  upgrade  the  UFPDB-2.0  model  for 
arbitrary  gate-dielectric  permittivity,  ed  = kde0,  where  the  kd  is  the 
dielectric  constant  and  e0  is  the  vacuum  permittivity.  In  the  upgraded 
code,  kde0  replaces  eox,  with  kd  being  a new  model  parameter. 

2.6  Conclusions 

Gate-to-body  tunneling  current  in  contemporary  MOSFETs  was 
physically  modeled  and  implemented  in  the  process-based  compact  model 
UFPDB  [Pel02],  which  is  unified  for  PD/SOI  as  well  as  bulk-silicon 
devices.  The  valence-band  electron  tunneling  (Jvbe)  to  the  n+-poly-Si  gate, 
predominant  for  inversion  conditions,  and  the  conduction-band  electron 
tunneling  (Jcbe)  to  the  substrate,  predominant  for  depletion/ 
accumulation  conditions,  were  modeled  by  modifying  the  classical 
independent-electron  formalism.  For  Jvbe>  several  scaling  effects  (i.e., 
quantization  effect  in  the  channel,  exchange  energy  of  inversion  electrons, 
and  Fermi-Dirac  statistics  in  the  high-doped  gate),  which  have  been 
overlooked  in  contemporary  modeling  of  gate  direct-tunneling  current, 
were  identified  and  accounted  for.  The  composite  IGB  model  showed  good 
agreement  with  measured  data  for  gate-oxide  thicknesses  varying  down  to 


1.65nm. 


CHAPTER  3 

ASSESSMENT  OF  THE  GATE-TO-BODY  DIRECT  TUNNELING- 
CURRENT  EFFECTS  IN  FLOATING-BODY  PARTILLY  DEPLETED  SOI 

CMOS  DEVICES  AND  CIRCUITS 

3.1  Introduction 

The  FB  effects  in  PD  SOI  MOSFETs,  and  the  resulting  hysteresis 
and  instability  in  dynamic  operations,  pose  major  challenges  for  very- 
large-scale-integration  (VLSI)  circuit  designs.  Further,  the  gate-to-body 
tunneling  current  modulates  the  body  voltage,  which  is  varied  by 
capacitive  coupling  in  transient  operation,  and  thus  can  influence  the  FB 
effects  [FunOO,  JosOl,  Chu02,  Poi02].  Such  influence  can  be  beneficial  or 
detrimental  for  the  dynamic  operations  of  FB  PD  SOI  circuits  and  must  be 
rigorously  examined  by  a reliable  physical  model.  The  dynamic  VGB  can  be 
positive  or  negative  in  dynamic  operation  of  FB  PD  SOI  devices  and 
circuits,  which  is  why  IGB  had  to  be  modeled  (in  chapter  2)  in  all  inversion, 
depletion,  and  accumulation  conditions. 

In  this  chapter,  we  use  an  IGB  model  in  UFPDB/Spice3  to  assess 
the  significance  of  IGB(t)  for  circuit  topologies  and  switching  patterns 
susceptible  to  FB  effects.  For  the  device/circuit  simulations,  the  UFPDB 
model  was  loosely  calibrated  to  a contemporary  scaled  device  with 
physical  gate  length  equal  to  of  60nm,  oxide  thickness  of  1.4nm,  and 
retrograded  channel  with  surface  doping  at  6.0xl017cm'3.  No  special 
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source/drain  engineering  for  carrier  recombination/generation  control  was 
assumed. 

3.2  Pass-Gate  Transistor  and  the  Transient  BJT  Current 
First,  we  consider  a pass-gate  transistor,  in  which  the  FB  effect 
tends  to  drive  the  transient  BJT  current  [Pel95],  possibly  resulting  in 
excessive  power  consumption,  degradation  of  noise  margin  and  stability, 
and  logic-state  error  in  pass-gate  and  other  dynamic-logic  circuits  [Chi96, 
Lu97,  Chu98],  We  assess  the  effect  of  reduced  transient  BJT  current  due 
to  the  gate-to-body  current  [Sin02].  Consider  the  FB  PD  SOI  nMOSFET 
with  source/drain  initially  biased  at  VGG  and  the  gate  grounded  in  a 
steady-state  condition  as  shown  in  Figure  3.1.  This  means  that  VGg  = - 
VDd>  with  VBS  = 0,  and  hence  holes  are  accumulated  in  floating  body. 
When  we  allow  for  gate-to-body  tunneling  current,  VBg  becomes  negative 
due  to  body  discharging  via  conduction-band  electron  tunneling  (IGbe); 
that  is,  the  body  voltage  is  now  less  than  VGB.  Now,  consider  the  source 
voltage  switching  to  ground  abruptly.  The  body  voltage  is  capacitively 
coupled  down  by  the  ratio  of  gate  capacitance  and  body-source  junction 
capacitance,  and  then  the  transient  VBg(t)  becomes  positive,  thereby 
driving  the  transient  BJT  current.  Since,  with  IGB,  the  body  voltage  is 
initially  less  than  VDD,  VBS(t)  is  decreased  relative  to  that  without  IGB, 
and  the  transient  BJT  current  is  reduced  as  shown  by  the  UFPDB/Spice3 
simulation  result  in  Figure  3.1,  by  about  50%  in  this  example  with  VDD  = 


1.0V. 
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Figure  3.1  UFPDB/Spice3-predicted  transient  BJT  current  in  the  PD/ 
SOI  pass  transistor  (nMOSFET:  Lgate  = 60nm,  tox  = 1.4nm) 
shown,  with  and  without  the  gate-to-body  tunneling  current 
(IGb  < 0)  accounted  for.  The  source-pulse  fall  time  is  0.2ns, 
with  VDD  = 1-0V. 
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Previous  work  concerning  PD  SOI  CMOS  pass-gate  circuits 
[Lu97,  Chu98]  has  shown  that  the  transient  BJT  current  can  be 
problematic  and  necessitate  modifications  in  circuit  design  to  control  its 
effects.  However,  as  the  CMOS  technology  is  scaled,  the  reduced 
tends  to  alleviate  the  BJT  effect  as  implied  by  the  discussion  concerning 
Figure  3.1.  Hence,  in  scaled  PD  SOI  CMOS  pass-gate  circuits,  we  do  not 
expect  the  gate-to-body  tunneling  current  to  provide  any  significant 
benefit,  even  though  it  can  reduce  the  transient  BJT  current  as  in  Figure 
3.1. 

3.3  Pass  Transistor-Based  Wide  Multiplexer 
To  check  the  Iqb  effect  on  pass  transistor-based  logic,  we 
examine  the  pass  transistor-based  n-to-1  multiplexer  circuit  in  Figure 
3.2(a),  where  n=32,  operating  at  VDD  = 1.0V.  We  consider  the  worst-case 
scenario  for  the  transient  BJT  effect  as  in  [Lu97].  All  inputs  are  assumed 
at  “1”  to  start  with.  The  control  signal  for  N32  and  P32  pair  is  “1”,  and 
hence  the  Vin  32  passes  “1”  to  node-1,  which  continues  to  be  “1”  afterwards. 
All  the  control  signals  (from  NT  and  PI  pair  to  N31  and  P31  pair)  are  “0”, 
and  the  (Vin  x - Vin  31)  are  then  pulled  down.  As  a result,  transient  BJT 
currents  flow  through  the  31  unselected  pairs  to  pull  down  node-1,  which 
is  supposed  to  be  at  “1”  and  is  held  by  the  single  selected  pair  (N32  and 
P32).  Unlike  the  previous  report  [Lu97],  which  shows  the  transient  BJT 
current  pulls  down  node-1  to  the  threshold  voltage  of  the  output  buffer  and 
can  lead  to  a wrong  logic  state,  the  UFPDB/Spice3  simulation  results  in 
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(a) 


(b) 

Figure  3.2  Pass  transistor-based  32-to-l  multiplexer,  (a)  Schematic  of 
the  circuit,  (b)  The  UFPDB/Spice3-predicted  effect  of  the 
gate-to-body  tunneling  current  on  the  node-1  voltage;  Lgate  = 
60nm,  tox  = 1.4nm,  and  Vdd  = 1-OV. 
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Figure  3.2(b)  show  that  the  node-1  voltage  is  not  pulled  down  enough  to 
upset  the  output  logic  state  even  if  we  consider  the  case  without  IGB.  This 
is  because  of  the  low  transient  B JT  current  due  to  low  VDD  as  discussed  in 
section  3.3.  Thus,  the  suppression  of  the  transient  BJT  effect  due  to  the 
Igb  is  n°t  relevant  in  this  case. 

3.4  Dynamic  Logic 
3.4.1  Dynamic  4-Way  OR  Circuit 

For  the  pass  transistor-based  wide  multiplexer  in  Figure  3.2, 
node-1,  which  is  pulled  down  by  transient  BJT  current,  is  held/restored  by 
the  single  selected  pass  gate.  The  effect  of  the  transient  BJT  current  is 
only  a very  small  dip  in  the  node-1,  and  thus  the  IGB  effect  is  not  relevant. 
For  dynamic  circuits  [Wes93],  the  effect  of  the  transient  BJT  current  can 
be  much  more  severe,  and  thus  the  IGB  effect  can  be  important  in  the 
circuit.  First,  consider  the  dynamic  four-way  OR  circuit  in  Figure  3.3(a), 
where  N1/N2/N3/N4  are  the  logic  transistors,  and  PO  and  NO  are  the 
clocked  precharge  and  evaluation  transistors,  respectively,  operating  at 
VDd  = 10V.  The  feedback  half-latch  pMOSFET,  which  holds  the  dynamic 
node-2  [Lu97],  is  eliminated  to  consider  the  worst-case  operation.  Assume 
in  a precharge  phase  that  the  input  to  N1  is  at  “1”  and  the  inputs  to  N2/ 
N3/N4  are  at  “0”  state.  The  dynamic  node-2  is  at  VBB,  and  the  common 
source  node-1  is  at  (VBB  - Vt).  The  input  to  N1  switches  at  t = 0.6ns  from 
“1”  to  “0”,  and  the  circuit  subsequently  evaluates  at  t = 1.1ns.  Transient 
BJT  currents  flow  through  the  off  devices  (N1/N2/N3/N4)  as  VI  drops  to 
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Figure  3.3  A dynamic  four-way  OR  circuit,  (a)  Schematic  of  the  circuit. 

(b)  The  UFPDB/Spice3-predicted  effect  of  the  gate-to-body 
tunneling  current  on  the  node-2  voltage.  The  width  of  each 
60nm  transistor  (with  tox  = 1.4nm)  is  set  to  maintain  the 
same  W/L  ratio  given  in  Figure  3(a)  of  [Lu97];  Vdd  = 1.0V. 
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“0”,  and  they  tend  to  pull  down  V2,  and  cause  an  erroneous  evaluation. 
The  UFPDB/Spice3  simulation  results  in  Figure  3.3(b)  show  that  the 
transient  BJT  currents  produce  a transient  voltage  dip  in  V2  to  0.85V 
without  IGB,  while  IGB  limits  this  dip  to  0.90V.  Because  of  the  low 
transient  BJT  current  due  to  low  VDD,  no  upset  occurs  (i.e.,  V2  is  not 
pulled  down  enough  to  cross  the  threshold  of  the  output  inverter)  in  either 
case,  but  we  note  that  IGB  does  result  in  less  dynamic  power  consumption 
and  more  noise  margin. 

3.4.2  Pseudo  Two-Phase  Dynamic  2-Way  NAND  Circuit 

The  transient  BJT  current  can  also  potentially  perturb  the  logic 
state  in  pseudo  two-phase  dynamic  logic  [Lu97]  as  shown  in  Figure  3.4(a). 
The  configuration  of  pseudo  two-phase  dynamic  logic  is  composed  of  two 
stages  of  dynamic  logic  whose  inputs  are  controlled  by  the  pass  gate  for 
each  stage.  The  pass  gates  are  controlled  individually  by  two  non- 
overlapping clocks  clkl  and  clk2.  While  the  first  stage  is  evaluating  by  “0” 
of  clkl  (“1”  of  clkl)  as  shown  in  Figure  3.4(a),  the  second  stage  is 
precharged  by  “1”  of  clk2  (“0”  of  clk2),  and  simultaneously  the  output  of 
the  first  stage  is  stored  in  the  gate  capacitances  of  the  second  stage  N3. 
The  subsequent  “0”  of  clk2  (“1”  of  clk2 ) forces  the  second  stage  to  be 
evaluated,  and  the  corresponding  output  is  latched  in  a succeeding  “1”  of 
clkl. 

The  time  period  when  clk2  is  at  “0”  is  considered  to  examine 
transient  BJT-current  and  corresponding  IGB  effects.  In  this  time  period, 
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Figure  3.4  A pseudo  two-phase  dynamic  2-way  circuit,  (a)  Schematic  of 
the  circuit,  (b)  The  timing  diagram;  Lgate  = 60nm,  Wgate  = 
5pm  (NO,  N3,  N4,  P0)  and  20pm  (Nl,  N2,  PI,  P2),  tox  = 1.4nm, 
and  Vj3D  = 1.0V. 
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the  pass  gate  to  the  second  stage  is  off  and  the  second  stage  is  evaluating. 
If  the  state  stored  in  the  gate  capacitance  of  N3  is  at  “1”  (node-3  is  at  “1”) 
and  the  first  stage  is  evaluated  to  be  “0”  (node-2  is  at  “0”)  at  the  falling 
edge  of  clkl,  transient  BJT  current  flows  through  the  nMOSFET  in  the  off- 
state  pass  gate,  and  the  current  discharges  the  input  node  (node-3)  from 
“1”  to  “0”  state.  The  transient  BJT  current  is  detrimental  to  the  circuit 
operation  via  the  discharge  of  the  input  node  voltage  for  dynamic  NAND 
logic.  A possible  example  is  illustrated  in  the  timing  diagram  in  Figure 
3.4(b).  During  the  evaluation  phase  for  the  second  stage  ( clk2  is  at  “0”), 
both  inputs  (node-3  and  node-4)  are  at  “1”  to  yield  a “0”  state  for  the  output 
Vout.  If  the  node-3,  which  is  held  by  the  gate  capacitance  of  N3,  is 
discharged  completely  by  the  transient  BJT  current  at  t = Ins  as  in  Figure 
3.4(b),  and  the  “1”  signal  arrives  at  node-4  at  t = 2ns,  the  NAND  logic  can 
be  erroneous  as  the  output  Vout  remains  in  the  logically  wrong  “1”  state. 
The  UFPDB/Spice3  simulation  results  in  Figure  3.5  show  that  the  node-3 
voltage  is  brought  down  to  0.70V  by  the  transient  BJT  currents  without 
Iqb>  while  Iqb  limits  this  to  0.80V.  Because  of  the  low  transient  BJT 
current  due  to  low  Vqq,  no  upset  occurs  in  either  case,  but  we  note  that 
Iqb  does  result  in  less  dynamic  power  consumption  and  more  noise 
margin. 

3.4.3  Dynamic  Cascade  Voltage  Switch  Logic  Circuit 

Dynamic  cascade  voltage  switch  logic  (dynamic  CVSL)  was 
investigated  as  another  circuit  topology  susceptible  to  the  transient  BJT 


60 


Time  (sec) 


Figure  3.5  The  UFPDB/Spice3-predicted  effect  of  the  gate-to-body 
tunneling  current  on  the  node-3  voltage:  Lgate  = 60nm,  tox  = 
1.4nm,  and  Vdd  = 1-0V. 
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effect  [Chu97],  We  show  how  the  IGB  can  influence  the  transient  BJT 
effects  in  the  circuit  using  UFPDB/Spice3.  Figure  3.6  shows  the 
schematics  of  a 3-input  dynamic  CVSL  XOR  circuit,  where  P0/P1  and  NO 
are  the  clocked  precharge  and  evaluation  transistors  respectively.  Pll/ 
Nil  and  P12/N12  are  the  output  inverters.  The  feedback  half-latches  to 
hold  the  dynamic  node-1  and  node-2  to  improve  the  noise  margin  in  Fig.  2 
of  [Chu97]  are  eliminated  to  consider  worst-case  operation.  The  width  of 
each  60nm  transistor  is  set  to  maintain  the  same  W/L  ratio  given  in  Fig. 
2(b)  of  [Chu97]  as  in  Figure  3.6. 

In  the  precharge  phase  ( elk  is  at  “0”),  the  clocked  evaluation 
transistor  NO  is  off  and  the  precharge  transistors  P0/P1  are  on.  For  any 
input  condition,  all  common-source  nodes  in  all  cascade  levels  are  at  “1” 
because  every  inputs  are  differential.  Therefore,  the  logic  transistors  with 
“0”  input  are  susceptible  to  transient  BJT  current  because  their  drains 
and  sources  are  at  “1”,  and  thus  the  bodies  are  at  “1”  if  enough  time  is 
given  to  yield  their  steady-state  values  before  the  circuit  goes  to  the 
evaluation  mode. 

In  the  evaluation  phase  ( elk  is  at  “1”),  with  the  input  pattern  (A, 
B,  C)  = (0,  1,  1),  the  node-3  is  pulled  down  by  the  clocked  evaluation 
transistor  NO.  Node-5  is  pulled  down  by  N2,  which  is  turned  on  (A  = 1). 
Transient  BJT  current  flows  through  the  Nl,  and  node-4  is  brought  down 
(supposedly  remains  at  “1”)  as  shown  in  Figure  3.7.  For  the  upper  level  of 
transistors,  node-7  is  pulled  down  by  N4  (B  = 1)  and  transient  BJT  current 
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Figure  3.6  Schematic  of  3-input  dynamic  CVSL  XOR  circuit;  Lgate  = 
60nm,  Wgate  = 7.5pm  (NO),  4.0pm  (N1  - N10),  0.6pm  (Nil, 
N12),  5.0pm  (Pll,  P12),  and  1.2pm  (P0,  PI),  tox  = 1.4nm,  and 
Vdd  = 10V.  “H”s  and  “L”s  are  the  logic  state  for  each  node 
and  the  arrows  represent  the  transient  BJT  current  in  the 
evaluation  phase  when  the  input  condition  is  (A,  B,  C)  = (0, 
1,  1). 
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Figure  3.7  The  UFPDB/Spice3-predicted  effect  of  the  gate-to-body 
tunneling  current  for  dynamic  CVSL  XOR  circuit:  Lgate  = 
60nm,  tox  = 1.3nm,  and  = 1.0V. 
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flows  through  the  N6  to  pull-down  node-6,  which  supposedly  remains  at 
“1”  state.  Note  that  pulling  down  node-7  by  turning-on  N4  induces  another 
transient  BJT  current  through  N5,  which  enhances  the  transient  voltage 
dip  of  node-4.  The  logic  transistors  of  C - C level  are  operated  by  the  same 
way  as  above.  Node-1  is  pulled  down  by  N8  and  transient  BJT  current 
flows  through  the  N10  because  node-7  goes  to  “0”  state.  The  transient  BJT 
current  pulls  down  node-2,  which  supposedly  remains  at  “1”  as  shown  in 
Figure  3.7.  Note  again  that  pulling  down  node-1  by  turning-on  N8  induces 
another  transient  BJT  current  through  N9.  As  can  be  seen  in  Figure  3.6, 
all  the  transient  BJT  currents  contribute  to  pull  down  node-2  (supposedly 
remains  at  “1”).  Logic  error  will  not  occur  because  the  transient  BJT 
current  is  not  large  enough  to  pull  down  node-2  across  the  threshold 
voltage  of  the  output  inverter  (P12/N12).  As  can  be  clearly  seen  in  Figure 
3.7,  however,  IqB  does  result  in  less  dynamic  power  consumption  and  more 
noise  margin  via  reducing  transient  BJT  current. 

3.5  Sense  Amplifier 

Finally,  we  consider  a sense  amplifier,  in  which  the  floating-body 
effects  are  typically  severe  and  for  which,  as  we  show,  IGB  can  provide  a 
benefit  even  in  scaled  PD  SOI  CMOS.  Without  IGB,  physical  circuit 
simulations  described  in  [Fos98],  based  on  a 0.25pm  CMOS  technology, 
showed  that  dynamic  instabilities  (i.e.,  bit  errors)  result  from  threshold- 
voltage  imbalances  caused  by  hysteretic  dynamic  floating-body  charging, 
necessitating  body  ties.  Here,  we  examine  the  effect  of  the  gate-to-body 
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tunneling  current  on  the  dynamic  instabilities  in  the  same  sense-amplifier 
circuit,  but  based  on  the  contemporary  Lgate  = 60nm  technology  noted 
previously.  The  schematic  of  the  full  CMOS  sense-amplifier  circuit  [Fos98] 
is  shown  in  Figure  3.8.  The  width  of  each  transistor  is  set  proportionally 
to  the  corresponding  W/L  ratio  used  in  the  original  design  [Fos98].  The 
sense  amplifier  is  composed  of  two  coupled  nMOS/pMOS  pairs,  Nl/Pl  and 
N2/P2.  We  assume  that  the  storage  capacitors,  e.g.,  CS1  and  CS2,  are  15 
fF,  and  the  bitline  capacitance,  represented  by  CBl,  is  250  fF. 

The  normal  sensing  operation  of  the  amplifier  is  indicated  by  the 
sequential  pulses  of  VPRE,  V(WL),  and  VSE  as  shown  in  Figure  3.8;  we 
assume  10ns  precharge  and  sense  times  with  0.5ns  rise  and  fall  times.  The 
sense  amplifier  must  respond  to  the  small  differential  voltage  (AV) 
established  across  the  bitlines  by  activation  of  WL1  (while  WL2  remains 
off).  For  example,  when  reading  a “0”  on  CS1,  AV  < 0 must  result  in  N1  and 
P2  being  turned  on,  and  N2  and  PI  off,  which  means  V(BL)  and  V(BL)  go 
to  0 and  VEE,  respectively.  The  bit  is  thus  read,  and  the  referenced  cell  is 
simultaneously  refreshed.  With  regard  to  possible  floating  body-induced 
instabilities,  NT  and  N2,  which  drive  the  proper  states  of  the  coupled 
transistor  pairs  in  response  to  AV,  are  the  most  crucial  devices. 

The  UFPDB/Spice3  simulation  results  in  Figure  3.9(a)  show  the 
body  and  source  voltages  of  N1/N2  during  a sequential  sense  cycle  in  the 
amplifier,  after  the  circuit  was  held  for  a lengthy  time  with  V(BL)  = 0 and 
V(BL)  = VEE  (as  in  an  extended  read-  or  write-“0”).  In  this  state,  VEg(N2) 
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Figure  3.8  DRAM  sense-amplifier  circuit,  including  precharge  and 
enable  circuitry  and  showing  two  complementary  data- 
storage cells  on  the  bitlines  and  representative  pulse 
sequences  for  sensing  DRAM-cell  data:  Lgate  = 60nm  and  the 
width  of  each  transistor  is  set  proportionally  to  the 
corresponding  W/L  ratio  used  in  the  original  design  [Fos98]. 
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Figure  3.9  UFPDB/Spice3  simulation  results  for  the  PD/SOI  DRAM 
sense  amplifier  (Lgate  = 60nm,  tox  = 1.4nm).  (a)  The  predicted 
NT  and  N2  body-voltage  transients  with  the  source  voltage, 
(b)  The  bit-line  voltage  with  and  without  gate-to-body 
tunneling  current  accounted  for.  For  Vqb  = 1.3V,  as  used 
here,  and  higher,  Iqb  prevents  the  erroneous  sensing  of  a “1”; 
for  lower  Vdq,  it  does  not. 
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= Vj)D  and  Vds(N1)  = 0,  and  hence  VBS(N2)  > VBS(N1)  as  defined  by  the 
the  body  nodal  equation  (IR(VBs)  = MvBd)  + !gb(vGb))  applied  to  both 
devices.  The  gate-to-body  current  (IGB)  of  N2,  which  is  conduction-band 
electron  tunneling  current  for  the  depletion  case  (Icbe  in  (2.26)),  is 
negligible  due  to  low  negative  VGB  (=  -VBS(N2)).  However,  IGB  of  Nl, 
which  is  valence-band  electron  tunneling  current  for  the  inversion  case 
(IvBE  in  (2.20)),  is  higher  due  to  high  VGg  (=  VBB),  and  can  significantly 
reduce  the  steady-state  and  subsequent  dynamic  differences  between 
Vbs(N2)  and  VBS(N1)  as  shown  in  Figure  3.9(a).  Hence,  for  high  enough 
VDD,  as  our  simulation  results  in  Figure  3.9(b)  show,  IGB  can  actually 
prevent  the  bit  error  resulting  from  the  floating-body  instability  noted  in 
[Fos98]. 

The  results  in  Figure  3.4  are  explained  as  follows.  In  the 
precharge  cycle,  the  drain  node  and  source  node  of  Nl  are  brought  up  from 
zero  to  VDD/2  and  VBB/2-Vt(t),  respectively,  since  Nl  is  ON.  The  drain- 
body-source  capacitive  coupling  in  Nl  causes  the  body  voltage  to  increase 
by  almost  VBG/2,  while  the  gate-to-body  capacitive  coupling  effect  in  Nl  is 
negligible  due  to  the  shielding  effect  of  the  inversion  layer.  However,  this 
is  not  the  case  for  N2,  where  there  is  no  inversion  layer.  The  drain  node  in 
N2  is  brought  down  to  VBG/2  and  the  gate  and  source  nodes  in  N2  are 
brought  up  to  VGD/2  and  VGG/2-Vt(t),  respectively.  The  net  capacitive 
coupling  effect  results  in  boosting  VB(N2)  up  to  0.82V,  irrespective  of  IGB. 
However,  IGB  in  Nl  makes  VBg(Nl)  higher  and  thus  decreases  the 
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imbalance  of  Vgg  in  N1  and  N2.  The  subsequent  dynamic  Vbs(N2)  > 
Vbs(NI)  imbalance  is  maintained  when  the  amplifier  is  enabled  at  t = 
50ns;  Vds(NO)  drops  to  0,  thereby  increasing  both  Vbs(N2)  and  Vbs(NI) 
via  the  gate-to-body-source  capacitive  coupling,  and  the  resulting  Vgs(t)’s 
define  unbalanced  Vt(t)’s.  This  dynamic  Vt(N2)  < Vt(Nl)  imbalance 
without  Iqb  is  large  enough  to  prevent  a valid  read-“0”  operation,  while 
the  threshold  voltage  imbalance  with  Iqb  is  not  enough  to  upset  the 
operation,  as  can  be  seen  in  Figure  3.9(b).  The  presence  of  Iqb  makes  the 
coupled  pairs  in  the  sense  amplifier  flip  properly  and  prevents  the 
erroneous  sensing  of  a “1”  state. 

All  of  the  simulations  discussed  in  this  section  did  not  account  for 
any  source/drain  engineering  to  enhance  carrier  recombination  and  effect 
good  control  of  the  floating-body  effects.  In  fact,  as  implied  in  Figure  3.9(a) 
by  the  high  DC  value  of  Vbs(N2)  in  the  sense  amp,  such  control  will  be 
necessary  in  scaled  PD/SOI  to  avoid  high  off-state  currents  without 
excessive  increase  in  threshold  voltage  [Pel02].  Additional  UFPDB/Spice3 
simulation  done  with  source/drain-junction  recombination  current 
increased  (by  more  than  an  order  of  magnitude)  show  that  all  the 
detrimental  floating-body  effects  are  well  controlled.  Hence,  we  conclude 
that,  in  general,  Iqb  should  be  insignificant  in  dynamic  operation  of 
optimally  scaled  floating-body  PD  SOI  CMOS  circuits. 
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3.6  Conclusions 

Using  the  Iqb  model  in  UFPDB/Spice3,  we  have  shown  how  the 
tunneling  current  tends  to  suppress  dynamic  floating-body  effects  in  some 
PD/SOI  circuits,  providing  benefits  that  follow  from  the  complex 
interactions  among  the  bias  conditions,  circuit  topologies,  and  switching 
patterns.  However,  the  significance  of  these  benefits  of  Iqb,  as  well  as  its 
detrimental  effects  [FunOO,  JosOl,  Chu02,  Poi02],  tends  to  be  diminished 
in  dynamic  operation  of  scaled  PD  SOI  CMOS  circuits  because  the  supply 
voltage  will  be  reduced  and  the  scaled  SOI  technology  will  necessarily  be 
optimized  to  suppress  the  detrimental  DC  floating-body  effects  on  off-state 
current. 


CHAPTER  4 

PARASITIC  BIPOLAR  JUNCTION  TRANSISTOR  IN  NONCLASSICAL 
ULTRA-THIN-BODY  TRANSISTORS:  MODELING,  VERIFICATION, 

AND  APPLICATION 

4.1  Introduction 

The  parasitic  BJT  effects  of  the  FB  SOI  MOSFET  have 
constituted  a critical  issue  for  a long  time  because  they  induce  a 
premature  breakdown  [You88,  Cho91],  latch  [Che88,  Cho91],  transient 
BJT  current  [Pel95],  and  also  amplify  GIDL  current  [Che92].  While 
impact  ionization-induced  phenomena  (premature  breakdown  and  latch) 
have  become  mitigated  due  to  scaled  operation  voltage,  the  transient  BJT 
current  and  the  amplified  GIDL  current  are  still  major  concerns  for 
nonclassical  UTB  transistors  regarding  control  of  off-state  leakage 
current  in  view  of  dynamic  and  static  power  consumption  [Sak03,  Fos99, 
Sch04] . 

However,  the  parasitic  BJT  model  [Kri96]  of  UFDG  [Fos03a], 
which  is  our  process/physics-based  compact  model  for  UTB  devices  and 
circuits,  is  inadequate  to  describe  the  transient  BJT  action  and  the 
amplification  of  GIDL  current  in  nonclassical  UTB  transistors  because  it 
does  not  correctly  account  for  accumulation  in  undoped  UTBs.  In  this 
chapter,  the  previous  BJT  model  is  modified  and  upgraded  to  correctly 
describe  majority-carrier  concentration  in  the  body,  which  governs  the 
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parasitic  BJT  gain.  Further,  circuit  application  with  the  upgraded  UFDG 
in  Spice3  [Fos03a]  is  exemplified  by  simulations  of  a nonclassical  dynamic 
logic  circuit,  which  show  the  parasitic  BJT  effect  on  circuit  operation  and 
stress  the  need  to  account  for  it  in  IC  design  at  future  CMOS  technology 
nodes. 

4,2  Effects  of  Device  Structure  on  Gate-Induced  Drain  Leakage  Current 

There  are  several  components  of  off-state  leakage  current  such 
as  gate-to-body/channel  tunneling,  band-to-band  tunneling  between  the 
body  and  drain,  direct  quantum  tunneling  from  source  to  drain,  drain- 
induced-barrier-lowering  (DIBL)-enhanced  subthreshold  conduction, 
GIDL  current,  and  gate  direct  tunneling  between  the  gate  and  source/ 
drain  extension  area  in  the  gate  overlap  region.  Among  these  leakage 
components,  GIDL  current  can  be  the  dominant  one  in  nonclassical  UTB 
transistors  [Sch04],  including  the  FDFET  and  the  double-gate  MOSFET 
(DGFET),  for  low-power  applications  because  it  can  be  enhanced  by 
increased  gate  work  function  needed  for  threshold-voltage  control,  and  it 
can  be  amplified  by  the  inherent  parasitic  BJT. 

For  FDFETs  and  symmetrical  DGFETs,  midgap  metal  is  needed 
for  gate  material,  instead  of  heavily  doped  polysilicon,  for  threshold- 
voltage  control,  while,  in  asymmetrical  DGFETs,  counter-type  polysilicon 
is  needed  for  the  back-gate  material.  Therefore,  the  gate  work  function  of 
nonclassical  UTB  transistors  will  be  larger  than  that  of  conventional 
MOSFETs,  and  thus  the  flat-band  voltage  of  the  gate-drain  MOS  structure 
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will  be  higher,  which  implies  higher  GIDL  current  due  to  higher  electric 
field  at  the  gate-drain  surface  interface.  With  the  assumption  that  the 
surface  potential  is  pinned  near  -Eg/q  + VBD  when  tunneling  occurs,  the 
normal  surface  field  is  expressed  as  [KimOl] 


Esd^ 


V VD  + - V 

VGD~VFB  + ~ V 


BD 


3 1 


oxf 


(4.1) 


The  higher  flat-band  voltage  of  the  nonclassical  UTB  transistors  makes 
the  surface  field  higher  (more  negative),  and  then  the  GIDL  current  is 
enhanced  via  the  following  conventional  model  [Fos99]: 


W = WA^|<-£.D)ex  (4.2) 

where  A and  B are  functions  of  the  bandgap  and  the  carrier  effective  mass, 
and  WALf  defines  the  effective  area  over  which  GIDL  occurs. 

Figure  4.1  shows  the  dependence  of  GIDL  current,  derived  from 
(4.1)  and  (4.2),  on  the  work-function  difference  between  the  gate  and 
drain-extension  region.  The  exponential  dependence  of  the  GIDL  current 
on  EsD  makes  the  strong  dependence  of  the  current  on  gate-oxide 
thickness  (toxf)  as  well  as  work-function  difference.  Note  that  the  work- 
function  differences  when  midgap  metal  and  counter-type  polysilicon  are 
used  as  gate  material  induce  huge  GIDL  currents.  Figure  4.2  shows  the 
effect  of  operation  voltage  (VDD  = -VGD)  on  GIDL  current  for  different  gate 
materials.  The  weak  dependence  of  the  current  on  VGd  suggests  that  even 
if  the  voltage  is  scaled  down,  GIDL  current  would  still  be  a dominant 
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Figure  4.1  The  predicted  dependence  of  GIDL  current  on  work  function 
difference  between  gate  and  drain  extension  region  in 

nMOSFET:  VGD  = -1.2V,  A = 2 x 109  A/cm2,  B = 4.5  x 109  V/ 
m,  and  ALf=10nm. 
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The  predicted  GIDL  current  vs.  operation  voltage:  EOT  = 
1.5nm,  A = 2 x 109  A/cm2,  B = 4.5  x 109  V/m,  ALf  = 10nm,  and 
VGD  = -VDD- 


Figure  4.2 
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component  of  off-state  leakage  current  in  nonclassical  UTB  transistors. 
Recently,  a simulation-based  study  showed  that  a gate-underlap  structure 
would  be  unavoidable  and  needed  for  nonclassical  UTB  transistors 
[Fos03b].  However,  even  for  the  gate-underlap  structure,  GIDL  current, 
which  can  be  induced  by  a gate  fringing  field,  can  still  limit  the  off- 
currents  for  low-standby  power  devices  [Sch04].  Therefore,  reliable 
physical  modeling  of  the  parasitic  BJT  is  needed  to  predict  the 
amplification  of  GIDL  current  in  nonclassical  UTB  transistors  employing 
metal  gates  and  undoped  bodies. 

4.3  Modeling  of  the  Parasitic  BJT 
From  the  condition  that  the  majority  carriers  are  approximately 
in  detailed  balance,  low-injection  parasitic-BJT  (base-transport)  current 
density  along  the  channel  (y->  direction  in  a generic  UTB  nMOSFET,  as 
illustrated  in  Figure  4.3,  is  expressed  (in  the  forward-active  mode)  as 


where  the  suffix  nL  stands  for  electron  and  low  injection,  respectively;  VT 
= kT/q  is  the  thermal  voltage.  Note  that  VBS,  the  body-source  voltage,  is 
defined  in  UFDG  via  the  body  nodal  equation  involving  recombination- 
generation  and  charging  currents.  (More  generally,  the  BJT  current 


[Kri96] 


(4.3) 
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Figure  4.3  Schematic  illustration  of  UTB  nMOSFET,  including  the 
parasitic  n+p'n+  BJT. 
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includes  an  exp(VBD/VT)  term  [Kri96],  as  does  our  subsequent  modeling  of 
it.) 

The  denominator  of  (4.3)  can  be  approximated  under  low- 
injection  conditions  to  be  [Kri96] 


f ' ff  p(x,  y)dy  = p(x)Leff 


(4.4) 


where  p(x)  is  the  gate-controlled  hole  density  in  the  UTB  at  depth  x where 
the  potential  is  minimum  along  y-direction;  Leff  is  the  effective  channel 
length  assumed  to  be  the  BJT  (quasi-neutral  in  y)  base  width.  The 
approximation  in  (4.4)  introduces  some  error,  which  can  be  tuned  out  via 
a model  parameter,  but  most  importantly  conveys  the  bias  dependence  of 
the  integrand.  Using  the  Boltzmann  approximation  for  carrier 
distribution,  the  hole  density  in  the  body  can  be  expressed  by 


p(x)  = Nb  exp 


ty(x)  -^BS 
VZ 


(4.5) 


where  is  the  body  doping  density  and  V|/(x)  is  the  potential  at  a point  x 
in  the  body,  referenced  to  that  of  a hypothetical,  unbiased  neutral  body. 
Note  that  the  front  surface  potential  (\|/sf)  is  assumed  in  the  previous 
model  [Kri96]  to  be  VBS  in  accumulation,  and  thus  the  hole  density  in 
accumulation  is  just  the  body  doping  density.  This  assumption  is  good  for 
a high-doped  body  (as  in  PD/SOI  MOSFETs  in  chapter  3).  However,  for 
undoped-UTB  transistors,  the  hole  density  in  accumulation  can  exceed  N^, 
and  thus  the  previous  BJT  model  largely  overestimates  the  parasitic  BJT 
current,  as  in  [Fos99].  Therefore,  the  potential  distribution  in  the  body  is 
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needed  for  accumulation,  as  well  as  weak  and  strong  inversion,  to 
calculate  p(x)  and  JnL  in  (4.3). 

For  inversion,  we  use  the  solution  of  potential  in  the  weak- 
inversion  model  [Yeh95]  of  UFDG  to  calculate  the  hole  concentration  in 
the  body,  and  also  use  the  same  model  for  strong  inversion,  but  pinning 
the  potential  to  (<\>q  (=  0.4V)  + ()>p)  [Tri04]  using  an  appropriate  smoothing 
function.  The  loose  validity  of  the  potential  model  in  strong  inversion  is 
mitigated  by  unnecessity  of  the  BJT  model  in  strong  inversion  where  the 
BJT  effects  are  insignificant  for  scaled  CMOS  technology  due  to  scaled 
V^P).  The  potential  model  in  strong  inversion  is  only  needed  for  numerical 
stability. 

In  order  to  calculate  the  potential  distribution  in  the  body  for 
accumulation,  we  relate  the  hole  concentration  to  the  accumulation  charge 
in  the  body  as 


tSi 


p(x)dx  = Qaccf 

(4.6) 

q(t%(x)dx  = Qaccb 

(4.7) 

2 


where  QaCcf  an<3  Qaccb  are  components  of  the  integrated  accumulation 
charge  in  the  body,  which  is  partitioned  at  the  center  as  indicated.  We 
assume  Qaccf  and  Qaccb  can  be  approximated  as  the  respective  gate 
charges: 

Qaccf  ~ -Coxf^GfS  ~ VbS  ~ ^ MS  + <t>C  “ 


(4.8) 
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Qaccb-~^'oxb^GbS~^BS~^>  AfS  + <|>c  - 0F)  . (4.9) 

Analogous  to  the  strong-inversion  assumption,  the  surface  potentials  (vysf 
in  (4.8)  and  \j/sb  in  (4.9))  are  estimated  to  be  pinned  at  (Vgg  -<| )q  + <J>p)  for 
strong  accumulation.  These  estimations  are  not  critical  because  the 
significant  dependences  on  the  actual  \)/sf  and  vj/gb  will  be  conveyed  by  the 
left-hand  sides  of  (4.6)  and  (4.7).  The  front  and  back  accumulation  charges 
in  (4.8)  and  (4.9)  are  reduced  to  zero  exponentially  with  increasing  gate 
bias  as  in  UFDG: 


Qaccf  = CoxfCln\  1 + exP 


( V GfS  ~ ^ MS  ~ V BS  + “ 0/iO 


C 


(4.10) 


Qaccb  = CoxbCln\  1 + eXP 


(^G6S  MS-  VbS  + _ 

C 


(4.11) 


where  C is  constant  which  defines  the  stiffness  of  the  smoothed  charge 
transition  between  accumulation  and  inversion. 

To  get  an  analytical  expression  of  (4.6)  and  (4.7),  first  for 
symmetrical  DGFETs.  the  potential  distribution  in  the  body  is 
approximated  as  linear  in  x between  the  surfaces  and  the  center  as  in 
Figure  4.4: 


V(*)  = 


NV-yc 


X 


for  0 < x < tSi/ 2 


(4.12) 


V|/(X) 


= Vc- 


yc-ySfr 

tSl/2  X 


for  tSi/2  < x < tSi . 


(4.13) 
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Figure  4.4  Potential,  \j/(x),  distribution  in  the  vertical  direction  and  the 
linear  approximation  for  symmetrical  double  gate  MOSFET. 
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Substituting  (4.5)  with  (4.12)  and  (4.13)  into  (4.6)  and  (4.7)  and  carrying 
out  the  integration  yields 


exp  - 


2 (lNbtSi 


Vc-vbs\  ( Vsf~VBS 


V , 


-exp  - 


_ vc  ^ 

vj  [ vj 


= Q 


accf 


(4.14) 


exp 


2 ^NbtSi 


Vsb  - VBS 

Vrr 


exp 


Vc-V 


BS 


Vr, 


Ysb) 

'VTJ 


Vc 


= Q 


accb 


(4.15) 


The  potential  \|/c  at  the  body  center  is  estimated  by  adding  aV-p,  which  is 
determined  by  numerical  simulation  and  adjusted  by  a model  parameter 
a,  to  the  average  value  of  the  surface  potentials: 


Vc 


2 


+ aVT  . 


(4.16) 


Using  (4.16)  and  equating  \j/sf  to  for  the  symmetrical  DG  device  then 
yields 


\\isf  - -VTln 


2a  Q 


accf 


.qNbtSi(l-e~a)_ 


+ vbs  = Vsb’ 


now  with 


Vc  = Vsf  + aVT- 


(4.17) 


(4.18) 


Note  that  (4.17)  and  (4.18)  analytically  relate  \\f  to  VQfg  = Vg^g. 

For  FDFETs  with  UTBs,  the  results  of  2-D  numerical  device 
simulations  with  MEDICI  [Med03]  in  Figure  4.5  suggest  that  \|/c  can  be 
assumed  to  equal  to  xj/gb,  for  varying  tgj,  at  least  for  long  Leff.  The  FD 


Potential  (V) 
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Vsb 


Figure  4.5  MEDICI-simulated  potential  distribution  (center  of  channel: 
y = 50nm)  of  FD  nMOSFET  with  varying  tgj  along  vertical 
direction  (x-direction);  Lmet  = lOOnm,  tox  = l.lnm,  tgox  = 
200nm,  midgap  gate,  NB  = 1.0xl0locm'3,  VGfg  = -0.5V,  VGbS 
= 0.0V,  and  VDS  = 1.2V. 
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nMOSFET  simulated  has  a long  channel  (Lmet  = lOOnm),  with  tgox  = 
200nm,  tox  = l.lnm,  and  undoped  body  (N^  = lxl010cm"3).  The  simulation 
results  shown  are  the  potential  \(/(x)  at  the  center  of  channel,  y = Leff/2. 
Carrying  out  the  integration  in  (4.7)  with  this  assumption  gives 


\\isb  - -Vrln 


'ZQaccb'- 

-<lNbtSi- 


+ VBS  = VC- 


(4.19) 


Now,  to  get  an  analytical  expression  for  \|/sf  for  the  FDFET,  the 
potential  distribution  is  approximated  as  constant  in  x between  the  front 
surface  and  the  center  as  illustrated  in  Figure  4.6.  This  approximation 
ultimately  yields  a high  \j/sf,  and  the  BJT  current  is  overpredicted  at  the 
front  surface.  However,  the  parasitic  BJT  current  is  predominant  between 
the  center  of  the  body  and  the  back  surface  due  to  relatively  low  p(x)  there. 
Thus,  the  error  due  to  the  approximation  in  Figure  4.6  is  negligible  for  the 
common  FDFET  structure,  which  has  a thick  buried  oxide.  Carrying  out 
the  integration  in  (4.6)  then  gives 


Vtf  = 


-V  T\n 


~2Qqcc/~ 

-Q^bhi- 


+ V 


BS‘ 


(4.20) 


Because  of  the  uncertainty  of  the  charge  model  in  (4.6)  - (4.9)  for  the 
FDFET,  a model  tuning  parameter  (a)  is  used  to  adjust  (4.19): 

¥sfc'(=  Vc')  = Vs&(=  Vc)~aVT  . (4.21) 

Hence,  \\i  is  analytically  related  to  Vqjs  and  VobS  f°r  th-e  FDFET. 
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Actual  potential 


Figure  4.6  Potential,  y(x),  distribution  in  the  vertical  direction  and  the 
constant  approximations  for  FDFET 
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Now,  for  both  the  DGFET  and  the  FDFET,  the  potential  in  the 
body  for  all  bias  regions  is  obtained  by  linking  the  potential  model  for 
accumulation  to  that  for  inversion  by  an  appropriate  smoothing  function: 

ln{  1 + exp[C(vj/ws  - V|/a)] } 


Vf  = Va  + 


(4.22) 


where  V|/a  represents  \|/sf,  \|/sb,  or  \| /c  for  accumulation  calculated  by  (4.16)  - 
(4.21),  and  \|/ws  represents  each  potential  in  weak  inversion,  which  is 
derived  from  the  weak-inversion  solution  [Yeh95]  of  UFDG.  Using  the 
solution  of  each  potential  (v|/t  in  (4.22))  in  the  body,  the  low-injection  BJT 
current  is  completely  described  by  (4.3)  - (4.5). 

The  high-injection  BJT  current  density  is  written  from  standard 
analysis  [Mul85], 


nH  ~ 


-2,D„n;exp(^ 

Leff 


(4.23) 


and  the  net  current  density  for  arbitrary  injection  level  is  given  as  [Kri96] 


T , \ ^nL^x)^nH 


(4.24) 


'H 


Substituting  (4.3),  (4.4),  and  (4.23)  in  (4.24)  yields 


J nTW  - 


o n 2 (VBS 
-elqDnni  expl — 


Jeff 


[2pix)+nieMwz  1 


V 


BS 


(4.25) 


The  total  BJT  current  is  calculated  by  integrating  (4.25)  between  x = 0 and 


x = tSi: 
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1 bjt  ~ -wf  J nj<{x)dx  (4.26) 

where  W is  the  device  width.  To  evaluate  (4.26),  a linear  variation  of 
potential  in  the  body  is  assumed  again  as  shown  in  Figure  4.4.  Using  (4.5), 
(4.12),  and  (4.13),  and  carrying  out  the  integration  in  (4.26)  with  (4.25) 
yields 


I 


BJT 


= w 


QDnn f 1 
Leff  PH 


{F1  + F2}exp 


(4.27) 


with 


i T m3  -i)*ri 

af  |_(fl-l)*3j 


(4.28) 


f2 


l,  r(^2- i)^3i 
bf\_(t3  - 1)^2 J ’ 


(4.29) 


u f^BS)  Vsf-Vc  u Vc-Vsfo  4.1  1,0/  r\\i 

where  p„  = ^exp^J,  a,  = -y-j-  , bf  = -y-j-  , tl  = 1 + 2p(x  = 0)/ 
PH,  t2  = 1 + 2p(x  = tgj)/pH,  and  t3  = 1 + 2p(x  = tgj/2)/pji.  When  \)/sf  equals 


\|/c  or  \|/c  equals  V|/Sb,  (4.28)  and  (4.29)  are  invalid  and  induce  numerical 
overflow.  In  that  case,  Fj  and  F2  in  (4.28)  and  (4.29)  are  smoothly  changed 
to  tgi/ti  and  tg;/t2  independently.  Eqation  (4.27)  defines  the  BJT  current 
for  forward  mode.  In  the  reverse  mode,  the  current  is  expressed  as  in 
(4.27),  with  Vbd  replacing  Vgg.  The  composite  parasitic  BJT  current  is 
then  given  by  the  summation  of  the  forward-mode  and  reverse-mode 


components. 
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For  the  transient  charge  dynamics,  the  bipolar  charges 
associated  with  each  of  the  terminals  must  be  properly  accounted  for  as  in 
[Kri96].  Of  the  two  main  components  of  stored  bipolar  charge  (minority 
carriers  stored  in  the  body  and  in  the  source),  the  previous  model  of 
minority  carriers  injected  in  the  body  should  be  modified  because  it 
depends  on  the  majority-carrier  distribution  in  the  body,  which  is  modified 
by  the  new  modeling  for  accumulation.  The  minority-carrier  charge 
injected  in  the  (p-type)  body  is  given  by 

Qn  = -W q^1  ^ff  \n(x,  y)dxdy  . (4.30) 

Assuming  a linear  distribution  of  electrons  along  the  channel/body  (base), 
for  all  injection  levels,  the  excess  electron  density  can  be  expressed  as 


An(x,  y)  = An(x)f  1 - ] . (4.31) 

^ Leff' 

Substituting  (4.31)  in  (4.30)  gives 

Qn^-WqI^(S,An(x)dx.  (4.32) 

2 Jq 


The  composite  An(x)  is  given  by  a weighted  sum  analogous  to  (4.24): 


An(x) = 


A nL(x)nH 
A nL(x)  + nH 


(4.33) 


where 


A nL(x) 


71: 


P(x). 


exp 


(4.34) 
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and  njj  = njexpCVgs  / 2Vr).  Using  (4.33)  and  (4.34)  in  (4.32)  and  carrying 
out  the  integration  yields  the  total  charge  stored  in  the  body  as 


WqLeff 


“■l°gJ 

af 


s2-fl 

sl-fl 


(4.35) 


where  fl  = exp(VBs  / V-p),  f2  = p(x  = 0)  exp(Vgg  / 2V-p)  / nj,  f3  = p(x  = t\J2) 
exp(Vgs  / 2Vt)  / nj,  SI  = fl  + £2,  S2  = fl  + f2  exp(aftb),  S3  = fl  + f3,  and  S4 
= f 1 + f3  exp(bft]3).  For  the  reverse  mode,  the  stored  charge  is  expressed  as 
in  (4.35)  with  Vbd  replacing  Vbs  in  Ph>  f2,  and  f3. 

4.4  Model  Implementation/Verification 

The  upgraded  model  for  the  parasitic  BJT  in  UTB  devices  was 
implemented  in  the  UFDG  model  in  Spice3.  An  existing  model  parameter, 
FVBJT,  which  was  previously  used  to  empirically  account  for  breakdown 
in  strong  inversion,  was  removed  since  the  parasitic  BJT  is  no  longer 
important  in  strong  inversion  due  to  scaled  operation  voltage.  A new 
model  parameter  FABJT  (a  in  (4.18)  for  the  symmetrical  DGFET  and  in 
(4.21)  for  the  FDFET)  was  added  for  tuning  the  BJT  gain  in  the 
accumulation  bias  condition.  Note  that  the  new  BJT  model  is  physically 
valid  only  for  symmetrical  DG  and  FD  MOSFETs.  For  asymmetrical  DG 
devices,  it  is  largely  empirical. 

To  predict  the  parasitic  BJT  gain  of  nonclassical  UTB  MOSFETs 
reliably,  the  calibration  of  the  smoothing-function  constant  C of  (4.22)  is 
critical  because  it  governs  the  body  potential  for  weak  accumulation, 
which  is  a common  bias  condition  related  to  GIDL  current  and  the 
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transient  BJT  current  of  a pass  gate.  The  smoothing-function  constant  is 
calibrated  via  numerical  simulation  of  the  transient  BJT  current  using 
MEDICI.  Table  4.1  gives  structural  information  for  the  simulated 


Table  4.1:  Device  structure  information 


Device  parameter 

DGFET 

FDFET 

Gate  length  (Lgate) 

48.0  [nm] 

48.0  [nm] 

Metallurgical  gate  length  (Lmet) 

28.0  [nm] 

28.0  [nm] 

Front  oxide  thickness  (toxf) 

2.0  [nm] 

1.1  [nm] 

Back  oxide  thickness  (t^) 

2.0  [nm] 

200.0  [nm] 

Silicon  film  thickness  (tgj) 

8.0  - 14.0  [nm] 

5.0  - 8.0  [nm] 

Body  doping  density  (N^) 

Intrinsic 

Intrinsic 

Source/Drain  doping  density  (Ngp) 

1020  [cm'3] 

1020  [cm'3] 

Source/Drain  doping  gradient 

1.0  [nm/dec.] 

1.0  [nm/dec.] 

Gate  work  function  (Oq) 

Midgap  + 0.1  [V] 

Midgap  + 0.1  [V] 

undoped-UTB  devices.  Lmet  (=  Leff)  was  defined  where  the  source/drain 
doping  density  drops  to  1015  cm'3.  For  the  transient  BJT-current 
simulation,  the  default  Auger  recombination  and  band-gap  narrowing 
models  in  MEDICI  were  used.  The  UFDG  model  parameters  related  to 
device  structure  were  directly  adapted  from  the  MEDICI-deflned 
structure.  (Note  that  quantization  effects  are  negligible  in  these  undoped 
devices  for  weak-inversion/accumulation  conditions.) 
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Figure  4.7  shows  UFDG-predicted  transient  BJT  current 
characteristics  of  the  symmetrical  DGFET  with  different  UTB  thickness, 
compared  with  MEDICI  predictions  for  the  pass-transistor  circuit.  Note 
that  the  dependence  of  parasitic  BJT  gain  on  tSi,  reflected  by  the  transient 
current  (which  was  almost  20  times  larger  in  the  previous  model  for  the 
same  device  structure)  is  modeled  well.  The  threshold  voltage  of  the  device 
with  thicker  tgi  is  lowered  by  SCEs,  and  hence  the  hole  density  in  the 
body,  which  defines  the  transient  BJT  current,  is  lowered.  Thus,  the 
symmetrical  DGFET  having  thicker  tgi  shows  higher  BJT  gain,  and  the 
upgraded  BJT  model  predicts  it  well. 

The  same  model  calibrated  to  the  DGFET  is  applied  to  the 
FDFET  with  device  parameters  in  Table  4.1  to  further  verify  the  new  BJT 
model.  The  model-predicted  transient  BJT  current  is  in  good  agreement 
with  the  results  of  MEDICI  simulations  as  shown  in  Figure  4.8.  The 
smoothing-function  constant,  C = 7.0  for  all  DG  and  FD  devices,  is  thus 
well  defined.  For  the  calibrations,  a in  (4.18)  for  the  symmetrical  DGFET 
and  in  (4.21)  for  the  FDFET  was  1.5  and  3.5,  respectively,  which  are 
reasonable  values.  The  physical  nature  of  the  new  BJT  model  is  thus 
confirmed. 

The  BJT  model  with  C evaluated  via  the  transient  BJT  current 
from  MEDICI  simulations,  is  next  applied  to  GIDL  current  amplification. 
The  measured  GIDL  data  of  UTB  FDFETs  [Kri02],  which  have  a midgap 
metal  gate  and  undoped  body,  are  used.  Calibrating  UFDG  [Fos03a]  to 
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Figure  4.7  UFDG-  and  MEDICI-predicted  transient  BJT  current  (see 
Figure  3.1)  for  symmetrical  DGFETs  having  several  UTB 
thicknesses.  The  supply  voltage  Vdd  was  1.2V  and  the  fall 
rate  of  the  source-voltage  pulse  was  1.2V/0.2ns.  Note  that  I^g 
is  normalized  with  respect  to  the  physical  width  of  the  DG 
devices. 
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Figure  4.8  UFDG-  and  MEDICI-predicted  transient  BJT  current  for 
FDFET  having  several  UTB  thicknesses.  The  supply  voltage 
VDD  was  1.2V  and  the  fall  rate  of  the  source-voltage  pulse 
was  1.2V/0.2ns. 
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measured  Ids'^GS  data,  and  tuning  a to  2.5,  we  model  the  GIDL  current, 
amplified  by  the  BJT,  to  be  quite  consistent  with  the  measured  data  for 
two  different  channel  lengths,  as  shown  in  Figure  4.9.  Hence,  the  BJT 
model  verification  is  solidified. 

The  parasitic  BJT  gain  implied  by  GIDL  current  also  depends  on 
silicon  film  thickness,  as  reflected  by  the  transient  BJT  current  in  Figures 
4.7  and  4.8.  The  validity  of  the  BJT  model  for  GIDL  current  versus  tgj  is 
examined  by  MEDICI  simulation.  Due  to  the  uncertainty  of  the  band-to- 
band  tunneling  model  in  MEDICI,  only  the  BJT  gain  (P),  calculated  by 
dividing  the  net  amplified  drain  current  by  the  GIDL  current  (which  is 
obtained  from  the  total  band-to-band  tunneling  current  in  the  output  file 
of  MEDICI),  is  compared  with  that  calculated  from  UFDG  simulation. 
Figure  4.10  shows  the  P(tgi)  results  of  UFDG  and  MEDICI  simulations  for 
FDFETs  and  symmetrical  DGFETs  with  Leff  = 28nm  at  two  different  gate 
biases.  The  UFDG  results  show  good  agreement  with  the  simulation 
results  of  MEDICI  at  weak  accumulation  (Vgs  = -0.15V)  for  both  devices, 
except  for  the  thickest  tgj.  This  discrepancy  is  related  to  the  low  threshold 
voltage  implied  by  worse  SCEs.  For  these  thick  devices,  the  channel 
current  is  comparable  to  the  GIDL-driven  BJT  current.  Thus,  the  parasitic 
BJT  gain  inferred  from  MEDICI  is  overpredicted  since  channel  current 
and  the  parasitic  BJT  current  cannot  be  distinguished  in  the  MEDICI- 
predicted  drain  current.  This  discrepancy  thus  does  not  necessarily  reflect 
model  error.  The  model  does,  however,  show  some  error  in  strong 
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Figure  4.9  Measured  [Kri02]  (points)  and  UFDG— predicted  current- 
voltage  characteristics  of  a UTB  FD  nMOSFET  showing  the 
parasitic  BJT  amplification,  (a)  Lgate  = 70nm,  toxf  = 2.1nm, 
tBox  = 200nm,  midgap  gate,  and  undoped  body,  (b)  Lgate  = 
80nm  and  the  other  parameters  are  same  as  in  (a).  The  tuned 
value  of  Bqjdl  is  3.8xl09V/m. 
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(b) 


Figure  4.10  UFDG-  and  MEDICI-predicted  parasitic  BJT  gain  for  GIDL 
current,  (a)  FDFET.  (b)  Symmetrical  DGFET. 
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accumulation  (Vqq  = -0.5V)  for  both  devices.  But,  the  strong-accumulation 
bias  condition  is  not  common  in  typical  operation  of  digital  logic  circuits, 
and  thus  this  deficiency  of  the  model  can  be  tolerable.  Hence,  overall,  the 
parasitic  BJT  model  is  effective  for  the  prediction  of  the  dependences  of 
BJT  gain  on  reasonable  bias  conditions  as  well  as  on  tgj  and  Leff. 

Note  that  the  BJT  gain  of  the  FDFET  is  much  higher  than  that 
of  the  DGFET  for  the  same  tgj  in  Figure  4.10.  For  the  FDFET,  the  low  hole 
concentration  at  the  back  surface  induces  large  parasitic  BJT  current. 
However,  tgj  of  the  FDFET  needed  for  SCE  control  is  thinner  than  that  of 
the  DGFET  for  the  same  Leff,  and  thus  the  BJT  gain  of  both  devices  with 
controlled  SCEs  will  not  differ  substantively. 

4.5  Model  Application 

To  examine  the  parasitic  BJT  effects  on  circuit  operation  for 
nonclassical  UTB  devices,  we  simulate  the  dynamic  4-way  OR  circuit  in 
Figure  3.3  of  chapter  3,  which  can  be  affected  by  the  transient  BJT 
current,  using  UFDG  in  Spice3.  The  schematic  of  the  circuit  is  redrawn  in 
Figure  4.11.  The  DGFET  for  the  simulation  has  Leff  = 28nm,  toxf  = toxb  = 
2.0nm,  intrinsic  body,  tgj  = 8nm,  and  midgap  gate.  For  the  FDFET,  Leff  = 
28nm,  toxf  = l.lnm,  toxb  = 200nm,  tgj  = 5nm,  intrinsic  body,  and  midgap 
gate  are  used.  The  width  of  each  transistor  is  set  proportionally  to  the 
corresponding  W/L  ratio  used  in  the  original  design  [Lu97].  As  assumed  in 
chapter  3,  the  dynamic  node-2  is  at  Vdd>  and  the  common  source  node-1  is 
at  VDD-Vt  in  precharge  phase.  When  the  circuit  goes  to  the  evaluation 
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Figure  4.11  Schematic  of  a dynamic  four-way  OR  circuit.  The  width  of 
each  28nm  transistor  is  set  to  maintain  the  same  W/L  ratio 
given  in  Figure  3(a)  of  [Lu97]. 
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mode  with  “0”  input  for  all  logic  devices  (N1/N2/N3/N4),  transient  BJT 
currents  flow  through  the  all  off  devices  (N1/N2/N3/N4)  as  VI  drops  to  “0”, 
and  they  pull  down  V2.  The  UFDG  simulation  results  for  both  DGFETs 
and  FDFETs  in  Figure  4.12  show  that  V2  is  pulled  down  by  the  transient 
BJT  current  enough  to  cross  the  threshold  of  the  output  inverter,  causing 
an  erroneous  evaluation.  On  the  contrary,  the  transient  BJT  currents  in 
the  FB  PD  SOI  circuit  were  not  enough  to  upset  valid  operation  due  to  low 
operation  voltage  in  chapter  3 (see  Figure  3.3(b)). 

Similar  to  the  results  in  Figure  4.7  and  4.8,  the  UFDG 
simulations  of  the  pass-gate  transistor  here  show  that  the  transient  BJT 
currents  in  the  DGFET  and  the  FDFET  are  comparable,  and  of  the  same 
order  of  magnitude  as  that  in  the  FB  PD  SOI  device  for  the  same  Vdd  (= 
1.0V).  Unlike  in  classical  CMOS,  however,  the  DGFET  and  FDFET  have 
very  low  intrinsic  gate  capacitance  for  low  VGS  (<  Vt)  [Fos02b,  Fos04].  For 
the  evaluation  mode  of  the  circuit,  the  node-2  voltage  is  mainly  held  by 
gate  capacitance  of  P5  and  N5.  Since  P5  is  in  the  off-state  (VGg  < Vt)  for 
the  evaluation  mode  and  the  output  inverter  is  skewed  to  speed-up  the 
output  pull-up  transition  (the  width  of  P5  is  8 times  larger  than  that  of 
N5),  the  total  capacitance  to  hold  the  node-2  voltage  is  much  smaller 
compared  to  that  of  the  FB  PD  SOI  circuit  (the  gate  capacitance  of  P5  is 
substantive  in  the  off-state  for  the  PD  SOI  device  [Fos02b]).  Thus,  the 
small  transient  BJT  current  can  pull  down  the  node-2  voltage  easily  for 
the  UTB  devices.  The  feedback  half-latch  is  needed  to  hold  the  dynamic 


100 


(a) 


(b) 

Figure  4.12  The  UFDG  predicted  effect  of  the  transient  BJT  current  on 
the  output  voltage:  Vdd  = 1-0V.  (a)  Symmetrical  DGFET.  (b) 
FDFET. 
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node-2  [Lu97],  and  the  size  of  the  half-latch  should  be  determined 
carefully  to  prevent  logic  upset  and  reduce  noise  margin  via  a reliable 
parasitic-BJT  model  for  the  UTB  devices. 

4.6  Conclusions 

The  effects  of  nonclassical  undoped  UTB  device  structures  on 
GIDL  current  have  been  examined.  GIDL  current  is  enhanced  by  the 
needed  higher  gate  work  function  of  the  UTB  devices  and  amplified  by  the 
parasitic  BJT.  To  predict  the  amplification  of  the  GIDL  current  as  well  as 
transient  BJT  current,  the  parasitic  BJT  model  for  nonclassical  UTB 
transistors  has  been  upgraded  and  implemented  in  the  UFDG  MOSFET 
model.  The  model  is  verified  with  the  results  of  MEDICI  simulations  and 
experimental  data  of  FDFETs  and  symmetrical  DGFETs,  and  shows  good 
agreement  for  various  device  dimensions  and  reasonable  bias  conditions. 
For  a dynamic  logic  circuit,  we  showed  that  the  parasitic  BJT  can  cause 
an  erroneous  evaluation  for  some  bias  conditions  and  switching  patterns 
due  to  inherent  low  gate  capacitance.  Thus,  a rigorous  simulation  is 
needed  to  ensure  right  operation  and  reduce  noise  margin  for  nonclassical 
UTB  transistor  circuits  via  a reliable  parasitic  BJT  model  for  the  UTB 


devices. 


CHAPTER  5 

TRIPLE-GATE  MOSFET:  INSIGHT,  FEASIBILITY,  AND  ANALYSIS 

5.1  Introduction 

Conventional  scaling  of  classical  bulk-silicon  and  PD  SOI  CMOS 
has  become  very  challenging  because  control  of  SCEs  for  acceptable  Ior/Ioff 
ratios  requires  precise  channel  doping  levels  and  gradients  [FraOl]  that 
are,  ultimately,  impossible  to  achieve.  Promising  alternatives  for 
continued  CMOS  scaling  are  the  nonclassical  devices  having  UTBs,  i.e., 
the  planar  FDFET  [Dor02]  and  the  quasi-planar  DG  FinFET  [Ked02], 
which  give  good  SCE  control  via  the  UTB  whether  it  is  doped  or  not.  The 
undoped  UTB  is  the  proper  choice  since  it  avoids  the  doping  challenges 
noted  for  the  classical  devices,  and  the  implied  threshold  voltage  (Vt) 
variations  [Xio03].  Further,  it  avoids  source/drain  junction-tunneling 
leakage  currents  and  can  perhaps  yield  high  carrier  mobilities  due  to 
negligible  impurity  scattering. 

However,  the  good  SCE  control  in  the  nonclassical  nanoscale 
devices  requires  that  the  UTB  thickness  (tgj)  be  considerably  thinner  than 
the  gate  length  (Lg),  as  well  as  being  uniform  across  the  device  and  circuit. 
For  the  FDFET  shown  in  Figure  5.1(a),  tgj,  which  is  the  height  (hgj)  of  the 
silicon  film,  should  be  scaled  down  to  ~Leff/5  [Tri03].  For  the  DG  FinFET 
shown  in  Figure  5.1(b),  tgj,  which  is  the  width  (wgj)  of  the  silicon  fin, 
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Figure  5.1  Schematic  illustrations  of  nonclassical  devices,  (a)  The 
FDFET.  (b)  The  DG  FinFET.  (c)  The  TGFET.  All  devices  are 
typically  fabricated  in  SOI  material,  with  the  buried  oxide 
(BOX)  providing  isolation. 
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should  be  scaled  down  to  ~Leff/2  [Fos03b],  or  about  the  minimum  feature 
size  in  conventional  CMOS  technology.  (Note  that  the  effective  channel 
length  Leff  of  nonclassical  devices  with  undoped  bodies  will  probably  be 
longer  than  Lg  [Fos03b].) 

The  TGFET  was  recently  proposed  [Doy03a,  Doy03b]  as  a means 
to  alleviate  the  stringent  thin-tgj  requirements  of  UTB  FDFETs  and  DG 
FinFETs,  enabling  “flexible  and  relaxed”  silicon-body  dimensions  as 
illustrated  in  Figure  5.1(c).  With  lower  hgj  and  wider  wgj,  the  planar-like 
TGFET  is  easier  to  fabricate  than  the  DG  FinFET,  and  the  experimental 
and  simulation  results  in  [Doy03a]  and  [Doy03b]  suggest  that  excellent 
short-channel  device  performance  can  be  achieved.  These  results, 
however,  were  obtained  for  TGFETs  with  heavily  doped  bodies,  which,  as 
noted  above,  are  not  desirable.  Further,  design  studies  of  TGFETs 
[Doy03b]  have  revealed  that  the  top  corners  of  the  silicon  body  (Figure  5.2) 
can,  with  strong  dependence  on  their  shape,  or  radius  of  curvature, 
significantly  affect  the  channel  current-voltage  characteristics  of  the  TG 
device.  Typically,  the  corner  components  of  current  reflect  a lower 
threshold  voltage  (Vt)  and  higher  off-state  current  (I0ff)  relative  to  the 
respective  properties  of  the  bulk  TG  device,  and  for  some  body  structures 
they  can  cause  a “hump”  in  the  subthreshold  log(Ii))-VG  characteristic 
[Doy03b].  Hence,  a serious  technological  issue  is  implied.  And,  the  layout 
area  of  the  TGFET  is  an  issue  for  nanoscale  CMOS  [Doy03b].  The 
feasibility  of  the  TGFET  is  thus  not  yet  clear. 
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Figure  5.2  2D  cross-sectional  view  A-A’  of  the  TGFET,  as  indicated  in 
Figure  5.1(c),  identifying  the  critical  corner  regions  of  the 
silicon  body/channel. 
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In  this  chapter,  we  examine  the  abnormal  corner  effects  in 
nanoscale  TGFET,  using  two-dimensional  (2-D)  numerical  device 
simulations  and  quasi-2-D  analysis  to  gain  physical  insights  on  how  the 
effects  could  possibly  be  suppressed,  irrespective  of  the  radius  of 
curvature,  which  is  difficult  to  control  technologically.  An  interesting  and 
effective  design  regarding  body/channel  doping  is  proposed  and  explained. 
Then,  we  examine  the  feasibility  of  the  TGFET  for  nanoscale  CMOS 
applications  with  regard  to  SCEs  and  layout  area.  We  check  the  SCEs  for 
varying  dimensions  of  the  body,  with  and  without  doping,  using  Davinci 
[Dav03],  a three-dimensional  (3-D)  numerical  device  simulator.  Then, 
with  the  SCEs  adequately  controlled  via  body  scaling  in  accord  with  the  3- 
D simulation  results,  we  analyze  the  gate  layout  area  of  the  integrated 
TGFET  needed  for  current  drive,  and  compare  it  with  those  of  the  FDFET 
and  DG  FinFET.  The  results  show  that  the  TGFET  is  not  feasible,  and 
that  the  DG  FinFET,  which  is  more  scalable  than  the  FDFET,  has  the 
most  potential  for  future  nanoscale  CMOS  applications. 

5.2  Suppression  of  Corner  Effects  in  TGFETs 
Scaled  TGFETs  are  3-D  devices  (Figure  5.1(c)),  and  hence  their 
complete  characterization  requires  3-D  analysis,  e.g.,  via  computationally 
intensive  3-D  numerical  simulation  as  in  [Doy03b].  However,  to  study  the 
basic  corner  effects,  we  can  analyze  the  2-D  cross  section  of  the  device, 
between  the  source  and  drain  as  shown  in  Figure  5.2,  which  is  tantamount 
to  assuming  long  gate  length  (Lg).  We  assume  a rectangular  SOI  body, 
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with  height  and  width  defined  as  hg^  and  wgg  respectively,  sitting  on  the 
underlying  buried  SiC>2  (BOX)/Si  substrate.  A reasonable  TGFET  design, 
with  a fully  depleted  body,  could  be  hgj  = wgj  = Lg  [Cha02,  Doy03b].  So,  we 
first  assume  a scaled  TG  nMOSFET  with  hgj  = wgj  = 30nm,  ignoring  the 
SCEs,  however,  as  noted  above.  We  let  the  gate  oxide  thickness  (tox  = EOT) 
be  l.lnm,  and  the  BOX  thickness  tgox  = 200nm.  Following  [Doy03b],  we 
further  assume  an  n+  polysilicon  gate,  and  heavy  doping  in  the  (p-type) 
body  for  Vt  control:  Ng  = 8.0xl018cm'3. 

Modeling  the  substrate  as  simply  a metal  contact,  which  is 
grounded,  we  use  MEDICI  [Med99]  to  solve  Poisson’s  equation  in  the 
defined  2-D  SOI  structure  (Figure  5.2)  for  various  gate  voltages,  Vq.  (For 
thick  tgQx,  effects  of  electric-field  fringing  and  induced  charge  in  the 
actual  silicon  substrate  are  small  and  can  be  ignored  here.)  Note  that  Vq 
equals  the  gate-source  voltage  (VQg)  of  the  TG  nMOSFET  when  the  source 
is  grounded.  A predicted,  1ow-Vq  2-D  inversion-electron  distribution  n(x, 
z)  in  the  body  cross-section  is  shown  in  Figure  5.3(a),  and  variations  of 
n(x=0,  z)  across  the  top  surface  (channel)  of  the  body  are  plotted,  versus 
Vq,  in  Figure  5.3(b).  Note,  for  the  lower  values  of  Vq  which  correspond  to 
weak  inversion,  the  dramatic  build-ups  of  the  electron  density  at  and  near 
the  top  corners  (z  = 0,  Wgj)  of  the  silicon  body.  They  reflect  the  reduced  Vt 
of  the  components  of  current  there  as  revealed  in  [Doy03b].  In  Figure 
5.3(b),  n(0,  z=15nm)  versus  Vq  reflects  a bulk  Vt  of  0.3-0.4V,  whereas  the 
corner  Vt  is  significantly  lower.  (Although  we  are  not  accounting  for 
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Figure  5.3  MEDICI-predicted  electron  concentration  of  a high-doped  TG 
nMOSFET.  (a)  2-D  electron  distribution  n(x,  z’),  where  z’  = 
z+100nm,  in  the  silicon-body/channel;  hgi  = wg^  = 30nm,  Vq 
= 0.3V.  (b)  Predicted  electron  density  n(x=0,  z)  across  the  top 
surface  of  the  cross  section  for  various  values  of  Vq.  Note,  for 

this  device  with  high  Ng  = 8.0  x 1018  cm'3,  the  extremely  high 
n in  the  corner  regions  for  low  Vq  (subthreshold  region).  This 
electron  buildup  for  high  Vq  is  much  less  significant. 
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quantization  (QM)  effects  [Ge02]  here,  our  results  concerning  the  corner 
effects  are  qualitatively  valid.  The  QM  solution,  which  requires  a 2-D 
Poisson-Schrodinger  solver  that  includes  anisotropic  effective  mass 
modeling,  would  show,  for  a specific  Vq,  smaller  n with  the  peak  densities 
displaced  from  the  surfaces,  but  still  with  relatively  high  densities  in  the 
corner  regions.) 

We  can  explain  the  inversion-charge  build-ups  at  the  corners  by 
considering  the  2-D  Poisson  equation  for  weak  inversion  (n  « Ng)  in  the 
p-type  silicon  body  (cross-section  in  Figure  5.2)  of  the  TG  nMOSFET: 

dEx  = qNB  dEz_  gNB{eff)  (5.1) 

dx  eSi  dz  eSi 

where  Ex  and  Ez  are  components  of  the  electric  field.  In  the  vicinity  of  the 
left-side  corner  (z  — » 0),  the  signs  of  both  partial  derivatives  in  (5.1)  are 
negative  for  a given  Vq  > VFB  (as  in  the  respective  1-D  versions  of  the 
equation).  Thus,  the  effect  of  the  left  gate  (i.e.,  Ez)  is  to  reduce  the 
magnitude  of  the  Ex  partial  derivative,  or  of  the  spatially  dependent 
effective  doping  density  (Ng(efg)  as  defined  in  (5.1).  This  reduction  of 
Ng(eff),  which  is  enhanced  for  z — > 0 (and  for  z — > Wgj),  defines  a reduction 
of  Vt  of  the  top  channel  near  the  corner(s).  (Note  that  if  wgj  were  large,  the 
1-D  Vt  of  the  top  channel  (with  Ez  = 0)  would  be  proportional  to  qNghgj 
[Lim83].  Also  note  that  our  explanation  of  diminishing  Vt  versus  z could 
be  given  analogously  in  terms  of  Vt  versus  x of  the  side  channels.)  The 
lower  Vt  in  the  corner  regions  is  consistent  with  thinking  of  these  regions 
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as  quasi- 1-D  symmetrical  DGFETs  having  an  effective  body  thickness 
(tsi(eff))  less  than  wgj,  where  Vt  depends  on  qNgt0Xtgi(eff)/2e0X  [Lim83]. 
Such  thinking  is  also  consistent  with  the  finding  in  [Doy03b]  that  the 
corner  devices  show  reduced  DIBL  since,  in  FDFET  and  DGFET,  the  SCEs 
tend  to  be  suppressed  as  tgj  is  scaled  down  [Tri03]. 

Because  of  technological  and  electrical  problems  associated  with 
controlled  heavy  doping  in  ultra-thin-silicon  bodies/channels,  viable 
nonclassical  devices  will,  more  than  likely,  be  undoped,  and  will  rely  on 
metal  gates  with  specific  work  functions  (Oq)  for  Vt  control  [Tri03]. 
Further,  we  surmise  that  low  Ng  in  (5.1)  could  influence  the  noted  corner 
effects.  So,  we  use  MEDICI  to  simulate  the  same  SOI  cross-section  of 
Figure  5.3,  but  with  Ng  = 0 (For  the  MEDICI  simulations,  we  set  Ng  = 
1.0xl010cnT3  to  avoid  numerical  instabilities.  Such  a low  Ng  (<  ~1015cm" 

Q 

) is  virtually  equivalent  to  Ng  = 0 with  regard  to  the  subthreshold 
characteristics  [Tri04].),  which,  for  a scaled  device  with  a small-volume 
body/channel,  will  be  the  actual  doping  density  in  the  undoped  device.  We 
assume  a midgap  metal  gate  (Oq  = 4.6V).  Predicted  n(x,  z)  versus  Vq  for 
this  device  is  shown  in  Figure  5.4.  The  bulk  Vt  is  a bit  higher  than  that  of 
Figure  5.3,  but  note  that  the  Iow-Vq  corner  effects  have  been  suppressed! 
For  high  Vq  (strong  inversion),  there  is  some  electron  build-up  in  the 
corner  regions,  but,  as  in  Figure  5.3,  it  is  relatively  insignificant,  in  accord 
with  results  in  [Doy03b].  In  this  case,  even  though  n > Ng,  (5.1)  with  Ng 
= 0 is  still  accurate  for  weak  inversion  because  the  inversion-charge 
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Figure  5.4  MEDICI-predicted  electron  concentration  of  an  undoped  TG 
nMOSFET.  (a)  2D  electron  distribution  n(x,  z’),  where  z’  = z 
+ lOOnm,  in  the  silicon-body/channel  cross  section  A-A’ 
(Figure  5.1(c));  hgj  = wgj  = 30nm,  Vq  = 0.4V.  (b)  Predicted 
electron  density  n(x=0,  z)  across  the  top  surface  of  the  cross 
section  for  various  values  of  Vq.  Note,  for  this  device  with  Ng 
= 0,  the  absence  of  any  n build-up  in  the  corner  regions  for 
low  Vq  (subthreshold  region).  The  electron  build-up  for  high 
Vq,  like  that  in  Figure  5.3,  is  not  significant. 
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(neglected  in  (5.1))  contribution  to  Vt  is  negligible.  Thus,  as  implied  by  the 
1-D  versions  of  (5.1),  both  the  Ex  and  Ez  partial  derivatives  are  relatively 
small,  and  hence  Ng(eff)  is  virtually  zero  and  there  is  no  significant 
reduction  in  Vt  as  the  corners  are  approached,  either  in  z or  in  x.  The 
nearly  uniform  n(x,  z)  in  Figure  5.4(a),  which  reflects  uniform  electric 
potential  in  the  silicon  body,  supports  this  explanation  of  no  corner  effects 
for  Ng  = 0.  And,  the  subthreshold  n(V(j)  variation  in  Figure  5.4(b)  reflects 
ideal  gate  swing,  S = 60mV/decade  at  T = 300K.  Additional  MEDICI 
simulations  of  TG  structures  with  hgj  and/or  wg;  increased  to  60nm  also 
showed  no  corner  effects  as  long  as  Ng  was  low.  We  note  that  such 
elimination  of  the  corner  effects  will  prevail  even  for  finite  Ng  defined  by 
natural  (unintentional)  doping  of  larger  bodies/channels. 

In  our  2-D  analysis  of  the  corner  effects  here,  we  have  ignored 
SCEs  (3-D  effects),  which  tend  to  reduce  Vt,  but  which  tend  to  be 
suppressed  by  scaling  the  silicon  body  dimensions  [Tri03].  So,  for  short  Lg, 
the  bulk  and  corner  Vt’s  will  be  lower,  but,  thinking  of  the  corner  devices 
as  symmetrical  DGFET  with  relatively  thin  tgj(eff)  as  previously  discussed, 
we  can  see  that  the  reduction  of  the  bulk  Vt  will  exceed  that  of  the  corner 
Vt.  Hence,  for  short  Lg,  the  corner  effects  in  the  high-Ng  device  will  be  less 
discernible  than  in  Figure  5.3.  However,  if  the  TGFET  is  properly  scaled 
to  control  SCEs,  this  amelioration  of  the  corner  effects  will  be  less 
significant.  The  need  for  ad  hoc  suppression  of  them  remains,  and  hence 
the  general  use  of  undoped  bodies  to  do  so  is  proposed.  This  proposal,  and 
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our  design  insights  regarding  scaled  TGFET,  are  supported  by  the  3-D 
simulations  of  corner  effects  in  scaled  FinFETs  and  shallow-trench- 
isolation  (STI)  bulk-silicon  MOSFETs  reported  in  [Bur03]. 

5.3  Short-Channel  Effects  of  Scaled  TGFETs 

With  the  insights  of  previous  section,  we  examine  the  feasibility 
of  the  TGFET  for  nanoscale  CMOS  applications  with  regard  to  SCEs  and 
layout  area.  We  check  the  SCEs  for  varying  dimensions  of  the  body,  with 
and  without  doping,  using  Davinci  [Dav03]. 

Following  [Doy03a],  we  first  assume  that  a reasonable  TGFET 
design  could,  with  reference  to  Figure  5.1(c),  be  hgj  = wgj  = Leff  (assumed 
to  equal  Lg).  For  Leff  = 28nm  devices,  we  let  the  gate-oxide  thickness  (tox 
= EOT)  be  l.lnm  and  the  buried-oxide  thickness  (tgox)  be  200nm.  For  the 
3-D  device  simulations,  we  assume  abrupt  source/drain  junctions  with 
lOnm  gate  overlaps  (Lg  = 48nm  with  Leff  = 28nm  being  the  metallurgical 
channel  length),  meaning  that  fringing  fields  outside  the  intrinsic  device 
are  negligible.  Figure  5.5  shows  Davinci-predicted  IdS'^GS  characteristics 
of  the  TG  nMOSFET,  with  both  doped  and  undoped  bodies.  For  the  doped 
device,  Ng  = 8.0xl018cm"3  and  the  gate  material  is  n+  polysilicon.  This 
device  exhibits  excellent  subthreshold  characteristics;  the  drain-induced 
barrier  lowering  (DIBL)  is  35mV/V  and  the  subthreshold  slope,  or  swing 
factor  (S)  is  75mV,  both  of  which  are  comparable  to  corresponding  results 
in  [Doy03a].  However,  as  we  have  previously  explained  in  previous 
section,  the  good  SCE  control  of  the  doped  TGFET  results  from  the 
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Figure  5.5  Davinci-predicted  current-voltage  characteristics  of  doped 
(Ng  = 8.0xl018cm'3  with  n+  polysilicon  gate)  and  undoped 
(midgap  gate)  n-channel  TGFETs;  tox  = l.lnm,  tgox  = 
200nm,  and  hgj  = wgj  = Leff  = 28nm. 
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predominant  subthreshold  current  flowing  in  the  corner  regions  of  the 
body,  which  can  be  thought  of  as  DG  devices  with  very  thin  effective  body 
thickness,  tSi(eff)  < wSi. 

To  examine  more  closely  the  doped-body  TGFET,  we  performed 
3-D  numerical  simulations  of  devices  with  varying  body  dimensions.  As 
shown  in  Figure  5.6,  the  predicted  subthreshold  characteristics  are 
virtually  independent  of  the  body  dimensions,  which  solidifies  our  insight 
that  the  subthreshold  characteristics  are  governed  by  the  corner  regions 
having  lower  Vt  than  that  of  bulk-body  device.  Note  in  Figure  5.6, 
however,  that  the  strong-inversion  channel  current  is  proportional  to  the 
effective  gate  width  (2hgj  + wgj)  of  each  device,  reflecting  three  surface 
channels  with  higher  conductance  than  the  corners. 

The  noted  insights  concerning  the  body-doped  TGFET  seem  to 
imply  an  optimal  device.  Corner  conduction  in  weak  inversion  yields  good 
control  of  SCEs  and  relatively  low  off-state  current  (I0ff),  whereas  the 
three  surface  channels  in  strong  inversion  imply  good  on-state  current 
(Ion),  as  shown  in  [Doy03a]  and  [Doy03b],  albeit  without  the  insights  noted 
here.  In  an  actual  TGFET,  however,  tgj(eff)  depends  on  the  finite  radius  of 
curvature  of  the  corners  [Doy03b],  which  could  be  difficult  to  control. 
Further,  controlled  doping  of  nanoscale  silicon  bodies  is  virtually 
impossible;  the  fluctuation  in  the  number  of  channel  dopants  and  their 
spacial  randomness  render  unacceptable  variations  in  Vt  [Xio03].  Hence, 
the  doped  TGFET  design  is  not  technologically  feasible. 
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Figure  5.6  Davinci-predicted  current-voltage  characteristics  of  doped- 
body  (Ng  = 8.0xl018cm'3  with  n+  polysilicon  gate)  n-channel 
TGFETs;  tox  = l.lnm,  tgox  = 200nm,  and  Leff  = 28nm.  (The 
numerical  simulation  of  the  largest-body  device  was  done 
using  larger  mesh-point  spacings,  which  could  yield  slightly 
erroneous  solutions  in  the  corner  regions  and  hence  underlie 
the  noticeably  lower  current  in  the  subthreshold  region.) 
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For  the  undoped  (NB  < ~lxl016cm‘3)  TGFET  in  Figure  5.5,  we 
assumed  a midgap  metal  gate  (Oq  = 4.6V)  for  threshold  control.  Note  that 
the  predicted  subthreshold  characteristics  are  much  worse  than  those  of 
the  doped  device  (and  are  prohibitive),  even  though  both  devices  have  the 
same  body  dimensions.  In  this  device,  the  corner  conduction  is  suppressed 
as  discussed  in  previous  section,  and  the  subthreshold  current  flows 
throughout  the  entire  body.  Obviously,  however,  the  body  is  too  large  to 
control  the  SCEs.  Nonetheless,  the  undoped  body  is  the  viable  design,  and 
hence  the  body  dimensions  needed  to  suppress  the  SCEs  must  be  checked. 

To  check  SCE  control  in  the  undoped-body  TGFET,  we  performed 
3-D  simulations  of  the  basic  structure  in  Figure  5.1,  still  assuming  Oq  = 
4.6V,  tox  = l.lnm,  tBOx  = 200nm,  and  abrupt  source/drain-body  junctions. 
Several  sets  of  body  dimensions  were  tried  in  the  simulations  to  obtain 
reasonable  subthreshold  characteristics,  i.e.,  DIBL  = lOOmV/V  and  S = 
80mV.  Required  combinations  of  hgj  and  wgj  for  the  TGFET  with  Leff  = 
28nm  are  shown  in  Figure  5.7;  in  fact,  the  dimensions  needed  are  about 
the  same  for  both  the  DIBL  and  S specifications.  Due  to  absence  of  the 
corner  effects  now,  the  body  dimensions  needed  to  suppress  the  SCEs  are 
much  smaller  than  those  intimated  in  [Doy03a],  e.g.,  hgj  = wgj  = Leff, 
indicated  by  the  ‘x’  in  Fig.  4,  which  yields  the  poor  subthreshold 
characteristics  shown  in  Figure  5.5.  Note  also  the  extreme  cases:  for  large 
hgj,  the  results  give  the  required  thin  wgj  for  the  DG  FinFET  (wgj(BG)  = 
0.53Leff)  to  get  the  specified  DIBL  and  S;  and  for  large  wSi,  they  give  the 
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Figure  5.7  Davinci-predicted  TGFET  body  dimensions  needed  (below 
and  left  of  the  curves)  to  suppress  SCEs  as  noted;  Leff  = 
28nm,  tox  = l.lnm,  VDD  = 1.0V.  For  the  extreme  cases, 
wSi(DG)  and  ^Si(FD)  represent  the  required  UTB  thinness  of 
the  DG  FinFET  and  the  planar  FDFET,  respectively,  to 
suppress  the  SCEs.  The  points  A and  B represent  the  body 
dimensions  for  hgj  = Leff  and  wgj  = Leff,  respectively,  showing 
that  TGFETs  have  negligible  scaling  advantage  over  the  DG 
FinFET  and  the  FDFET.  Point  C is  an  “optimal”  TGFET 
design,  having  wg;  and  hgj  equal  to  about  1.4  times  wgjQQ) 
and  hgj(pD),  respectively.  Note  that  the  design  with  hgj  = wgj 
= Leff  (point  x)  is  clearly  not  feasible. 
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required  thin  hSi  for  the  FDFET  (hSi(FD)  = 0.21Leff).  These  UTB 
requirements  are  consistent  with  results  in  [Tri03]  and  [Fos03b]. 

To  stress  the  much  more  stringent  body  scaling  needed  for  the 
nanoscale  TGFETs  relative  to  that  (hgj  ~ wgj  ~ Leff)  implied  in  [Doy03a] 
and  [Doy03b],  we  note  in  Figure  5.7  (point  A)  that  with  hgj  = Leff,  wgj  = 
wgi(DG)  is  required  to  suppress  the  SCEs.  This  means  that  the  TGFET  has 
no  scaling  advantage  over  the  DG  FinFET  at  all!  Clearly  then,  hSi  should 
be  made  as  high  as  possible  to  maximize  the  effective  gate  width,  creating 
in  essence  the  DG  FinFET  with  gate  width  ~2hgj.  Also,  we  note  in  Figure 
5.7  (point  B)  that  with  wSi  = Leff,  the  required  hSi  for  the  TGFET  is  less 
than  20%  larger  than  that  needed  for  the  FDFET  (hgi(FD)).  This  means 
that  the  TGFET  has  very  little  scaling  advantage  over  the  FDFET! 

Our  results  and  conclusions  regarding  the  TGFET  versus  the 
FDFET  are  inconsistent  with  recently  reported  experimental  results 
[Kri03],  which  suggest  that  TGFETs  with  wSi/Lg  >2.8  show  excellent 
subthreshold  characteristics  down  to  Lg  ~ 25nm,  even  for  hgj/Lg  ~ 0.4.  Our 
results  in  Figure  5.7  show  that  hSi/Lg  < 0.2  would  be  needed  to  get  the 
results  reported  in  [Kri03].  Whereas  there  could  be  uncertainty  in  Leffvs. 
Lg,  and/or  in  tgj  (=  hgj),  in  [Kri03],  the  authors  suggest  that  local  strain  in 
the  body  due  to  the  metal  (NiSi)  gate  and/or  the  SOI  mesa  isolation 
benefits  the  SCE  control  (as  well  as  carrier  mobility  and  Ion).  Strained 
TGFETs  might  then  be  feasible,  but  understanding  and  controlling  the 


strain  constitute  a formidable  task. 
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5.4  Layout  Area 

With  the  required  body  dimensions  for  SCE  control  given  in  the 
preceding  section,  we  now  examine  the  CMOS  layout-area  issue  of 
TGFETs,  and  compare  it  with  those  of  planar  FDFETs  and  DG  FinFETs. 
We  focus  on  gate  areas,  assuming  that  the  peripheral  areas  of  the  three 
devices  scale  proportionally.  Multi-finger  structures  as  illustrated  in 
Figure  5.8,  with  their  pertinent  dimensions,  will  generally  be  required. 
The  2-D  cross-sectional  view  A-A’  of  the  multi-gate/finger  transistor, 
indicated  in  Figure  5.8  (a),  is  shown  in  Figure  5.8  (b)  to  relate  the  effective 
gate  width  (Weff)  of  each  finger  (per  pitch  in  the  layout)  to  the  pitch  (P).  In 
order  to  calculate  the  gate  layout  areas  of  the  TGFET  and  the  DG  FinFET 
needed  to  achieve  the  same  total  drive  current,  we  assume  Ion  per  unit 
width  of  each  device  is  the  same,  and  equal  to  that  of  the  FDFET,  which 
will  be  used  as  a reference.  With  this  assumption,  gate  layout  areas  for 
multi-gate/finger  TGFETs  and  DG  FinFETs  are  expressed  as 


A-tgfet  ~ Lg 


- P 

■2hSi  + wSi  J 


(5.2) 


A FinFET  ~ Lg 


r W. 


&Si 


(5.3) 


where  Wg  is  the  total  gate  width  needed  (which  defines  the  number  of 
fingers).  Note  that  (5.2)  and  (5.3)  depend  on  Weff  for  each  device;  i.e.,  as 
indicated  in  Figure  5.8  (b),  Weff  = 2hgj  + wg;  for  the  TGFET  and  Weff  = 2hg; 
for  the  DG  FinFET.  For  the  planar  FDFET,  the  gate  layout  area  (Ajtd)  is 
simply  WgLg,  and  Weff  = P. 
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Figure  5.8  The  schematic  illustrations  of  the  multi-gate/finger 
transistor,  (a)  Top  view,  (b)  Cross-sectional  view  A-A’  as 
indicated  in  (a).  The  effective  gate  width  of  each  finger  (per 
pitch  in  the  layout)  is  noted  for  the  TGFET,  the  DG  FinFET, 
and  the  planar  FDFET  in  (b). 
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Based  on  the  simulation  results  in  Figure  5.7,  we  assume,  for 
SCE  control,  wgj/Leff  = 0.5  for  the  DG  FinFET  and  hgj/Leff  = 0.2  for  the 
FDFET.  The  purported  advantage  of  the  TGFET  is  the  relaxation  of  the 
silicon-body  dimensions  required,  relative  to  the  DG  FinFET  and  the 
FDFET.  However,  this  advantage  is  severely  restricted  due  to  SCEs  as 
illustrated  in  Figure  5.7.  In  order  to  relax  wgj  and  hgj  at  the  same  rate 
while  maintaining  the  same  SCEs,  the  body  dimensions  of  an  “optimal” 
TGFET  are  chosen  to  be  about  1.4  times  wgjQQ)  and  hgi(FD),  as  indicated 
by  point  C in  Figure  5.7.  Then,  using  (5.2)  and  (5.3),  we  show  in  Figure  5.9 
calculated  ratios  (Ra)  of  A^gfet  and  Afreet  t°  Apo  versus  Lg  (assumed  to 
equal  Leff)  at  different  technology  nodes,  with  all  body  dimensions  set  as 
noted  for  SCE  control.  For  the  calculations,  Lg  and  P were  obtained  from 
the  ITRS  [SemOl]  CMOS  projections  for  the  high-performance  (HP)  logic 
technology  and  the  low  standby-power  (LSTP)  technology.  We  assumed 
maximum  aspect  ratios  (Rf  = hgjAvgj)  of  5 for  the  DG  FinFET.  Since  Ra  of 
the  DG  FinFET  is  inversely  proportional  to  Rf,  it  would  increase  in  Figure 
5.9  by  a factor  of  5/3  ~ 5/4  for  smaller  Rf  = 3 ~ 4,  which  is  perhaps  more 
technologically  reasonable  for  the  scaled  device. 

As  evident  in  Figure  5.9,  Ra  of  the  TGFET  is  much  higher  than 
that  of  the  DG  FinFET,  which  can  be  <1  for  higher  Rf.  At  the  Lg  = 37nm 
node  of  the  HP  logic  technology,  the  TGFET  would  require  almost  four 
times  the  gate  area  required  for  the  DG  FinFET.  This  severe  layout-area 
disadvantage  prevails,  for  the  same  body-dimension  factor  of  1.4,  all  the 
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Lg  (nm) 


Figure  5.9  Calculated  layout-area  ratios  (Ra)  of  the  TGFET  with  wg}  = 
1.4wgj(QQ)  and  hg;  = 1.4hgj(pD)  and  the  DG  FinFET  with  Rf  = 
hgi/wgi  = 5 relative  to  the  planar  FDFET,  versus  gate  length; 
Lg  (=  Leff)  and  P for  the  calculations  were  obtained  from  the 
ITRS  projections  for  HP  and  LSTP  CMOS  technologies 
[SemOl]. 
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way  to  the  end  of  the  ITRS  where  Lg  ~ lOnm.  Note  also  in  Figure  5.9  that 
Ra  is  better  for  both  devices  in  the  LSTP  technology  because,  at  a given 
node  for  specified  P,  Lg  is  scaled  less  aggressively.  (For  a specific  Lg,  the 
HP  vs.  LSTP  differences  in  Figure  5.9  reflect  different  values  of  P.) 
However,  the  severe  TGFET  layout-area  disadvantage  remains.  Finally, 
we  note  that  all  the  ratios  plotted  in  Figure  5.9  would  be  reduced  by  a 
factor  of  two  if  the  multi-gate/finger  devices  were  defined  by  spacer 
lithography  [Doy03b],  thereby  substantially  enhancing  the  relative 
layout-area  advantage  of  DG  FinFETs. 

The  layout-area  comparisons  in  Figure  5.9  reflect  the  smaller 
Weff  for  the  TGFET  structure  used.  To  check  its  optimization,  the  3-D 
simulation  results  for  S = 80mV  in  the  TGFET  with  Leff  = 28nm  in  Figure 
5.7  are  replotted  in  Figure  5.10.  Lines  A and  B there  represent  Weff  = P for 
the  TGFET  and  DG  FinFET,  respectively,  as  for  the  planar  FDFET.  To  get 
Weff  equal  to  or  wider  than  P,  the  respective  body  dimensions  must  be 
above  the  lines  A and  B.  For  the  TGFET,  our  “optimal”  design  (point  C in 
Figure  5.7)  is  well  below  line  A.  We  could  redesign  it  as  indicated  in  Figure 
5.10  for  wider  Weff  and  smaller  Ra.  However,  it  is  evident  that  moving 
toward  line  A with  the  same  SCE  control  would  require  either  decreasing 
wgi  (path  ® in  Figure  5.10),  which  yields  a virtual  DG  FinFET  (Weff  of  the 
TGFET  is  wider  than  that  of  the  DG  FinFET  having  the  same  body 
dimensions  due  to  active  top  gate,  but  the  added  gate  width  is  relatively 
small  when  hg/wgj  >>  1.),  or  decreasing  hgi  (path  ® in  Figure  5.10),  which 


125 


Figure  5.10  Davinci-predicted  body  dimensions  of  the  TGFET  needed  for 
S = 80mV;  Leff  = 28nm,  tox  = l.lnm,  and  VDD  = 1.0V,  with 
lines  A and  B representing  Weff  = P for  the  TGFET  and  the 
DG  FinFET,  respectively.  The  arrows  indicate  possible 
design-optimization  paths,  but  lead  to  DG  FinFET  and 
planar  FDFET  structures. 
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yields  a virtual  FDFET.  The  TGFET  structure  is  inferior,  unless  it  is 
transformed  into  a DG  FinFET  or  a FDFET!  Note  also  in  Figure  5.10  that 
the  DG  FinFET  would  be  made  more  area-efficient  than  the  FDFET,  i.e., 
have  a wider  Weff,  if  hgj  could  be  above  line  B;  this  would  require  Rf  > 5. 

From  the  3-D  device  simulation  results  in  Figure  5.7,  we  have 
noted  the  requirement  of  body  thickness  for  SCE  control  to  be  wgj/Leff  - 
0.5  for  the  DG  FinFET  and  hSi/Leff  ~ 0.2  for  the  FDFET.  Hence,  if  we  set 
the  pragmatic  lower  limit  of  body  thickness  (wgj  or  hgj)  as  5nm  due  to 
manufacturing  burden  and  quantization  effects  from  the  structural 
confinement  [Tri03],  we  note  that  the  FDFET  cannot  be  scaled  down  below 
Lg  ~ 25nm  (where  we  are  still  assuming  Leff  ~ Lg.).  However,  the  DG 
FinFET  can  be  scaled  down  to  Lg  - lOnm.  Indeed,  this  superior  scalability 
of  the  undoped-body  DG  FinFET  seems  to  render  it  the  most  promising 
candidate  for  future  nanoscale  CMOS. 

5.5  Conclusions 

Via  2-D  numerical  device  simulations  and  quasi-2-D  analyses  of 
nanoscale  TGFETs,  we  have  shown  and  explained  that  the  reduced  Vt  of 
the  corner  regions  in  the  body/channel  can  be  eliminated  by  leaving  the 
body  undoped,  and  hence  relying  on  a metal  gate  with  proper  work 
function  for  Vt  control.  This  finding  adds  to  the  technological  and 
electrical  reasons  [Tri03]  for  proposing  the  use  of  undoped  bodies  in 
nonclassical  MOSFETs;  the  problem  of  controlling  the  shape  of  the  corners 
in  the  TGFET  [Doy03b],  which  is  probably  not  possible,  is  eliminated. 
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Focusing  on  viable  undoped  bodies,  we  examined  the  UTB 
dimensions  needed  to  control  SCEs  in  nanoscale  TGFETs  by  3-D 
numerical  simulations.  The  results  showed  that  much  more  stringent  body 
scaling  is  needed  for  undoped  TGFETs,  which  have  no  corner  effects, 
relative  to  the  doped  ones,  which  are  technologically  and  electrically 
infeasible.  When  the  undoped-body  dimensions  are  scaled  for  adequate 
SCE  control,  we  find  that  the  TGFET  suffers  from  a significant  layout- 
area  disadvantage  relative  to  the  DG  FinFET  and  the  planar  FDFET.  We 
hence  conclude  that  the  TGFET  structure,  without  some  benefit  from 
strain  [Kri03],  is  not  feasible  because,  when  adequately  scaled  for  SCE 
control,  its  gate  layout-area  efficiency  is  relatively  poor  and  can  be 
improved  only  by  evolving  it  into  a DG  FinFET  with  hSi/wSi  >>  1 or  a 
FDFET  with  hg/wgj  <<  1.  Further,  our  insights  imply  that  this  layout 
disadvantage  is  not  peculiar  only  to  the  TGFET  but  is  characteristic  of  any 
multi-gate  MOSFET  with  more  than  two  gates,  for  example,  the  gate-all- 
around  device  [Col90]  and  the  surrounding-gate  transistor  [Tak91].  We 
thus  argue  that  the  DG  FinFET,  which  is  more  scalable  than  the  FDFET, 
is  the  nonclassical  device  with  the  most  potential  for  future  nanoscale 
CMOS  applications. 


CHAPTER  6 

SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  WORK 

6.1  Summary 

In  this  dissertation,  analysis  and  modeling  of  parasitic  effects  in 
advanced  SOI  CMOS  technologies,  including  UTB  transistors  and 
TGFETs,  were  presented,  and  the  feasibility  of  the  TGFET  was  discussed. 
The  major  contributions  of  the  research  are  summarized  as  follows. 

In  chapter  2,  gate-to-body  tunneling  current  (IGB)  in 
contemporary  MOSFETs  was  physically  modeled  and  implemented  in  the 
process-based  compact  model  UFPDB.  The  valence-band  electron 
tunneling  (Jvbe)  to  the  n+-polysilicon  gate  of  the  nMOSFET,  predominant 
for  inversion  conditions,  and  the  conduction-band  electron  tunneling 
(Jcbe)  to  the  substrate,  predominant  for  depletion/accumulation 
conditions,  were  modeled  by  modifying  the  classical  independent-electron 
formalism.  For  Jvbe>  several  scaling  effects  (i.e.,  quantization  effect  in  the 
channel,  exchange  energy  of  inversion  electrons,  and  Fermi-Dirac 
statistics  in  the  high-doped  gate),  which  have  been  overlooked  in 
contemporary  modeling  of  gate  direct-tunneling  current,  were  identified 
and  accounted  for.  The  composite  Igg  model,  applicable  to  pMOSFET  as 
well,  showed  good  agreement  with  measured  data  for  gate-oxide 
thicknesses  varying  down  to  1.65nm. 
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In  chapter  3,  using  the  Iqb  model  in  UFPDB/Spice3,  we  showed 
how  the  tunneling  current  tends  to  suppress  dynamic  floating-body  effects 
in  some  PD/SOI  circuits,  providing  benefits  that  follow  from  the  complex 
interactions  among  the  bias  conditions,  circuit  topologies,  and  switching 
patterns.  However,  the  significance  of  these  benefits  of  Iqb>  as  well  as  its 
detrimental  effects  [FunOO,  JosOl,  Chu02,  Poi02],  tend  to  be  diminished  in 
dynamic  operation  of  scaled  PD/SOI  CMOS  circuits  because  the  supply 
voltage  will  be  reduced  and  the  scaled  SOI  technology  will  necessarily  be 
optimized  to  suppress  the  detrimental  DC  floating-body  effects  on  off-state 
current. 

In  chapter  4,  the  effects  of  nonclassical  undoped  UTB  device 
structure  on  GIDL  current  were  examined.  GIDL  current  is  enhanced  by 
the  needed  higher  gate  work  function  of  the  UTB  device  and  amplified  by 
parasitic  BJT.  To  predict  the  amplification  of  the  GIDL  current  as  well  as 
transient  BJT  current,  the  parasitic  BJT  model  for  nonclassical  UTB 
transistor  was  upgraded  and  implemented  in  the  UFDG  MOSFET  models. 
The  model  was  verified  with  the  results  of  MEDICI  simulations  and 
experimental  data  of  FDFETs  and  symmetrical  DGFETs,  showing  good 
agreement  for  various  device  dimensions  and  reasonable  bias  conditions. 
For  a dynamic  logic  circuit,  we  showed  that  the  parasitic  BJT  can  cause 
an  erroneous  evaluation  for  some  bias  conditions  and  switching  patterns 
due  to  inherent  low  gate  capacitance.  Thus,  a rigorous  simulation  is 
needed  to  ensure  right  operation  and  reduce  the  noise  margin  for 
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nonclassical  UTB  transistor  circuits  via  a reliable  parasitic  BJT  model 
dependent  on  the  UTB  device  structure. 

In  chapter  5,  we  showed  and  explained  that  the  reduced  Vt  of  the 
corner  regions  in  the  body/channel  of  triple-gate  FinFETs  (TGFETs)  can 
be  eliminated  by  leaving  the  body  undoped,  and  hence  relying  on  a metal 
gate  with  proper  work  function  for  Vt  control.  This  finding  adds  to  the 
technological  and  electrical  reasons  [Tri03]  for  proposing  the  use  of 
undoped  bodies  in  nonclassical  MOSFETs;  the  problem  of  controlling  the 
shape  of  the  corners  in  the  TGFET  [Doy03b],  which  is  probably  not 
possible,  was  eliminated.  Focusing  on  viable  undoped  bodies,  we  examined 
the  UTB  dimensions  needed  to  control  SCEs  in  nanoscale  TGFETs  by  3-D 
numerical  simulations.  The  results  showed  that  much  more  stringent  body 
scaling  is  needed  for  undoped  TGFETs,  which  have  no  corner  effects, 
relative  to  the  doped  ones,  which  are  technologically  and  electrically 
infeasible.  When  the  undoped-body  dimensions  are  scaled  for  adequate 
SCE  control,  we  found  that  the  TGFET  suffers  from  a significant  layout- 
area  disadvantage  relative  to  the  DG  FinFET  and  the  planar  FDFET.  We 
hence  conclude  that  the  TGFET  structure,  without  some  benefit  from 
strain  [Kri03],  is  not  feasible  because,  when  adequately  scaled  for  SCE 
control,  its  gate  layout-area  efficiency  is  relatively  poor  and  can  be 
improved  only  by  evolving  it  into  a DG  FinFET  with  hSi/wSi  » 1 or  a 
FDFET  with  hg/wgj  <<  1.  Further,  our  insights  imply  that  this  layout 
disadvantage  is  not  peculiar  only  to  the  TGFET  but  is  characteristic  of  any 
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multi-gate  MOSFET  with  more  than  two  gates,  for  example,  the  gate-all- 
around  device  [Col90]  and  the  surrounding-gate  transistor  [Tak91].  We 
thus  argue  that  the  DG  FinFET,  which  is  more  scalable  than  the  FDFET, 
is  the  nonclassical  device  with  the  most  potential  for  future  nanoscale 
CMOS  applications. 

6.2  Suggestions  for  Future  Work 
The  following  research  tasks  are  suggested  as  future  work 
regarding  advanced  SOI  devices. 

For  a high-K  dielectric,  the  gate-to-channel  tunneling  current 
would  still  be  considered  as  a major  leakage  current.  Modeling  and 
analysis  of  the  gate-to-channel  tunneling  current  are  suggested  for  UTB 
transistors,  with  (or  without)  high-K  dielectrics. 

Recently,  a study  based  on  simulation  shows  that  a gate- 
underlap  structure  would  be  unavoidable  and  needed  for  nonclassical  UTB 
transistors.  However,  even  for  the  gate-underlap  structure,  a numerical 
simulation-based  study  shows  GIDL  current,  which  can  be  induced  by  the 
gate  fringing  field,  can  still  limit  the  off-currents  for  low-standby  power 
devices.  Therefore,  the  modeling  of  GIDL  current  in  nonclassical  UTB 
transistors  employing  the  gate-underlap  structure  is  recommended. 

2-D  device  simulation  shows  that  unamplified  GIDL  current 
depends  on  UTB  thickness  in  symmetrical  DGFETs,  i.e.,  GIDL  current 
increases  for  thinner  bodies.  This  dependence  perhaps  results  from 
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volume  inversion  in  the  thin  Si  film  under  the  gate-drain  overlap. 
Modeling  of  GIDL  current  in  UTB  devices  is  thus  recommended. 


APPENDIX 

DERIVATION  OF  THE  INTERBAND  TUNNELING  CURRENT 

EQUATION 

Previous  studies  of  direct  tunneling  currents  in  MOS  devices 
have  focused  primarily  on  the  component  which  flows  between  the  gate 
and  the  channel,  i.e.,  intraband  (conduction  band-to-conduction  band) 
tunneling.  However,  modeling  of  the  tunneling  current  that  flows  between 
the  gate  and  the  body  (as  in  chapter  2),  i.e.,  interband  (valence  band-to- 
conduction  band)  tunneling,  has  not  been  adequately  documented.  In  the 
interband  tunneling,  the  proper  characterization  of  the  tunneling  current 
is  somewhat  different  from  the  commonly  used  direct  tunneling  current 
expression  [Sha99,  Dep95,  She72]  since  the  electron  density-of-states 
effective  mass  is  different  in  the  valence-  and  conduction-bands  and  the  Si 
bandgap  is  indirect,  which  implies  the  involvement  of  phonons  in  the 
tunneling  process  [Mol64],  However,  we  do  not  explicitly  account  for  the 
conservation  of  transverse  momentum  (via  phonons)  in  this  case;  in  fact, 
the  functional  dependence  of  tunneling  probability  on  electric  field  is 
almost  identical  [Mol64]  to  that  for  a direct  bandgap  semiconductor,  even 
though  the  needed  phonon  involvement  decreases  tunneling  probability 
dramatically  [Mol64].  With  this  insight,  we  develop  here  a general 
expression  for  both  interband  and  intraband  tunneling  current  by 
generalizing  the  classical  model  of  Harrison  [Har61]. 
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The  general  probability  per  unit  time  of  the  tunneling  transition 
of  an  electron  in  a state  a on  one  side  of  the  tunneling  region  to  a state  b 
on  the  other  side  was  given  as  [Har61] 


2k 


Pab  “ 1^|Mab|  Pb^a^1-^) 


(A.l) 


|Mab| : the  matrix  element  for  the  transition 
pb : the  density  of  states  at  b 

fa,  fb:  the  probabilities  of  occupation  of  the  states  a and  b,  respectively. 
The  same  formalism  applies  to  the  reverse  tunneling  process: 


Pb,  = f |Mba|2pafb(l-fa). 


(A.2) 


The  total  tunneling  current  is  derived  by  summing  over  all  states  a of 
fixed  kt  (transverse  momentum),  summing  over  kt,  multiplying  by  2 for 
spin,  and  multiplying  by  the  electronic  charge  q: 


Jab  = 2qXCPabPadE> 


(A.3) 


Jba  = 2qirPbaPbdEx. 


Hb 


The  net  current  is  then 


j=jab-Jba  = 2q(lJ>abPadEx-Xl>baPbdE3 


^tb 


(A.4) 


(A.5) 


By  replacing  the  sum  over  kt  by  an  integral  over  a constant 
energy  surface,  i.e., 
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I = L2J^12)tktdkt  = L2j^_dEt 

K V^r  J (2tc)2  J 2 TTh 


(A.6) 


which  follows  from 


E.  = 


, 2 
K kt 

2mt* 


for  2D  electrons,  we  have,  per  unit  area  (L2), 


= 2q^[mta*£‘,"“XPabP,dExdEt-mtb*f“'^pbapbdExdE; 
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Mab|  PbPa^a(1  _^b)^x^t 


-m 


tb 


.(A.7) 


With  Harrison’s  formalism  [Har61]  for  the  matrix  element  in  (A.7), 


|Mab|2paPb  = —2e~" 

4k 


(A. 8) 


f^ox 

where  rj  = 2J^  |kx|dx  (with  kx  being  the  momentum  perpendicular  to  the 


Si02-Si  interface), 


(A.7)  becomes 


-m 


tb 


'Jo*  "W  Jo  e ’'fbtl  - fa)dExdEt]  (A.9) 


which  is  a general  (interband  or  intraband)  tunneling  equation. 
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